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ABSTRACT 


When  a  semiconductor  device  fails,  a  complicated  set  of 
actions  takes  place  to  move  the  device  from  a  high  voltage, 
low  current  state  to  a  low  voltage,  high  current  state. 
Accompanying  these  changes  are  the  collapse  of  isotropic 
current  flow  to  a  small  column  or  filament  of  current  which 
rapidly  grows  and  heats  the  device  to  the  melting  point. 

These  actions  are  collectively  called  second  breakdown. 

This  work  reviews  past  results  to  show  that  second  breakdown 
is  the  formation  of  a  current  controlled  negative  resistance 
CCCNR5  regime  which  necessarily  forces  current  flow  to  form  a 
current  filament.  The  current  filament  is  modeled  as  a  solid 


state  plasma  column  undergoing  a  self-induced  magnetic  pinch 
as  a  result  of  the  CCNR .  The  dispersion  relation  is  derived 
to  first  order  to  show  that  the  pinching  leads  to  an  unstable 
equilibrium  that  could  lead  to  material  failure. 

The  minimum  set  of  parameters  necessary  to  the  formation 
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CHAPTER  I 


THE  PROBLEM  AND  ITS  HISTORICAL  BACKGROUND 

At  first  glance  mast  scientists  ujould  expect  that 
determining  the  conditions  necessary  to  cause  a  transistor  or 
diode  to  fail  luould  be  related  to  a  feiu  simple  factors. 

These  factors  ujould  most  certsinly  involve  the  current  and 
voltage.  Failure  mould  then  he  caused  by  literally  melting 
the  device  with  a  large  amount  of  current  or  forcing  the 
device  to  arc  over  from  a  high  voltage. 

Indeed,  there  is  very  little  proof,  considering  the 
material  these  devices  are  made  of,  that  the  above  explana¬ 
tion  of  device  failure  is  untrue.  Common  experience  shows 
that  when  a  diode  is  forced  to  pass  several  amperes  of 
current,  when  it  is  rated  for  only  a  fraction  of  an  ampere, 
something  is  going  to  melt.  Similar  situations  can  be 
constructed  for  Just  about  any  semiconductor  device  available 
today . 

In  the  early  1360 ’s  the  United  States  Air  Force  CUSAF) 
and  several  private  firms  attempted  to  determine  exactly  what 
conditions  of  current  and  voltage  caused  semiconductor 
devices  to  fail.  These  first  researchers  took  the  obvious 
path  already  described  above  and  found  evidence  that  devices 
failed  at  levels  very  different  than  expected.  In  addition, 
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the  devices  exhibited  phenomena  that  no  one  could  consist¬ 


ently  explain. 

This  mark  luili  attempt  to  explain  and  determine  the 
conditions  under  luhich  a  semiconductor  device  uill  exhibit 
signs  of  impending  failure  due  to  an  electrical  overstress. 
Device  failure  uill  be  linked  to  the  onset  of  negative 
resistance,  which  will  be  shown  to  force  current  in  the 
device  to  collapse.  Th.is  collapse  will  be  modeled  as  a  solid 
state  plasma  that  undergoes  a  self-induced  magnetic  pinch  to 
an  unstable  equilibrium  under  negative  resistance  conditions. 
The  resultant  instability  will  also  be  shown  to  be  capable  of 
inducing  rapid  device  heating. 

These  results  will  lead  to  the  quantitative  development 
of  a  minimum  pinch  criteria  that  will  be  compared  to  a  set  of 
measurements  on  a  transistor  that  is  forced  into  the  condi¬ 
tions  that  are  known  to  precede  device  failure.  The 
remainder  of  this  chapter  will  detail  the  background  of 
device  failure  research  and  how  experimental  evidence  led 
this  author  to  conclude  that  device  failure,  and  the  forma¬ 
tion  of  a  negative  resistance  regime  in  the  device,  were  one 
and  the  same. 

The  Phusics  of  Device  Failure:  The  Problem 

Historu  of  Research  Results 

The  first  published  reports  concerning  device  failure 
were  reported  for  diodes  by  Tauc  and  Abraham  Cll  in  1957. 


Failure  of  transistors  ujas  soon  reported  by  Thornton  and 
Simmons  CSD  in  1950.  These  reports  demonstrated  experimen¬ 
tally  that  the  failure  state  of  a  device  was  characterized  by 
a  transition  from  a  high  voltage,  iou  current  state,  to  a 
high  current,  louj  voltage  state.  Soon  after  this  transition 
the  device  failed  as  an  open  circuit  or  short.  Autopsies  of 
failed  devices  shouied  def.’nitive  signs  of  extensive  heating 
leading  to  melting  of  the  semiconductor  material . 

For  semiconductor  devices  made  of  silicon  or  germanium, 
this  implied  that  temperatures  were  approaching  the  melting 
point  of  the  materials  during  failure.  Hany  researchers 
concluded  that  such  temperatures  mere  a  result  of  ohmic 
heating  of  the  device  in  this  new  regime  Cof  high  current, 
low  voltage)  that  occurred  after  avalanche  breakdown.  Thus, 
many  researchers  referred  tc  this  regime,  which  occurred 
after  avalanche  or  first  breakdown,  as  second  breakdown. 

Continued  research  soon  broke  into  two  approaches  in  an 
attempt  to  determine  device  failure  levels  C33 .  One  group 
approached  failure  as  a  solely  thermal  event.  A  second  group 
sought  further  information  concerning  why  a  device  went  into 
this  new  regime  of  operation  called  second  breakdown.  In  the 
early  1960 ’s,  a  tremendous  number  of  articles  were  published 
concerning  device  failure,  of  which  a  partial  listing  is  set 
forth  in  reference  three. 
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Tuio  very  general  conclusions  came  forth  from  these  reports: 

1.  Attempts  to  model  device  failure  as  a  thermal 
runaway  event  C4I1  had  very  limited  success  and 
suffered  from  extreme  problems  of  repeatability 
:5: . 

E.  Devices  that  entered  second  breakdown  suddenly 
exhibited  the  collapse  of  isotropic  current  flow  to 
a  small  current  filament  or  column  IBD .  This 
filament  soon  engulfed  the  device  volume,  melting 
the  device  as  it  grew  L7J . 

Such  results  were  quite  unexpected  and  continued  to  fuel 
further  research  along  both  lines  of  interest. 

Researchers  who  followed  conclusion  one  above  made 
numerous  attempts  to  model  device  failure  by  utilizing 
extensive  thermodynamic  approaches  coupled  with  device 
geometry.  These  initial  ideas  became  the  basis  for  further 
research  and  were  developed  into  models  by  Ulunsch  and  Bell 
CBl  in  1966  and  Tasca  CS]  in  1970.  These  models  generally 
utilized  a  thermal  diffusion  model  incorporating  the  increas¬ 
ing  carrier  concentration  of  a  semiconductor  with  increasing 
temperature  driven  by  an  external  current.  The  feedback 
mechanism  provided  by  a  semiconductor’s  increasing  carrier 
concentration,  with  temperature  which,  in  turn,  increases  the 
current,  became  the  prime  physical  explanation  of  device 
failure  in  these  models.  The  two  models  developed  by  UJunsch 
and  Tasca  could  not  incorporate  all  aspects  of  device  design 
necessary  for  an  accurate  model,  so  a  correction  factor  was 
added  to  each  model  by  matching  its  results  empirically  to 
experiment . 
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Initially,  results  at  applying  tnese  models  to  simple 
diodes,  and  transistors  treated  as  diodes,  looked  very 
promising  (see  Figures  1  and  2) .  The  only  information 
required  to  determine  a  device’s  failure  level  consisted  of 
the  material  parameters  and  the  appropriate  device  area. 
Further  application  of  these  tico  models  to  other  diodes  and 
transistors  shoued  very  limited  results,  and  more  important¬ 
ly,  the  correction  factor  incorporated  in  the  models  uiouid 
have  to  be  changed  In  an  almnsc  random  fashion  to  bring 
predictions  closer  to  experimental  results  CIOD.  Perusal  of 
these  studies  shou  that  some  results  uere  dissimilar  due  to 
differences  in  the  definition  of  what  constituted  device 
failuire.  Variation  of  what  constituted  the  onset  of  failure 
also  presented  problems  in  interpreting  experimental  ap¬ 
proaches  used  by  different  researchers. 

The  community  of  researchers  trying  to  refine  thermal 
models  tried  to  clear  up  confusion  by  utilizing  the  onset  cf 
second  breakdown  as  the  criterion  for  device  failure. 

Several  researchers  then  tried  tc  alter  the  old  thermal 
models  to  explain  second  break.down  as  an  electronic  and 
thermal  event.  These  electrer hermai  mcdeis  C11,1SI'  were,  by 
themselves,  achievements  in  modeling  device  dynamics  which 
occurred  on  totally  different  time  scales  Cthermal  events 
occur  in  microseconds,  electrical  events  occur  in 
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Figure  1.  Experimenta.!  data  points  for  the 
IN459  diode.  Slope  indicates  pouier  is  inversely 
proportional  to  the  square  r-oot  of  the  time  to 
failure  as  dictated  by  the  Uunsch-Bell  model. 

The  diode  is  pulsed  in  reverse  bias  and  failure 
corresponds  to  the  device  assuming  a  resistor 
characteristic.  Only  seven  devices  were  tested. 
From  reference  8. 
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Figure  2.  A  composite  of  data  points  For 
eight  transistors.  All  failure  uas  by  pulsing 
the  emitter  base  Junction  in  reverse  bias  with 
the  collector  open.  Lines  indicate  the  theoret¬ 
ical  failure  levels  predicted  by  the  UJunsch-Bell 
model.  Failure  coresponds  to  the  device 
assuming  a  resistor  characteristic.  From 
reference  0. 
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times  on  the  order  of  nanoseconds).  The  general  results  of 
these  models  did  show  better  agreement  with  experiment,  but, 
they  were  limited  by  the  extensive  information  required  of  a 
device  to  determine  failure.  In  simple  terms,  the  model  was 
as  accurate  as  the  device  information  was.  Since  the 
electrothermal  models  required  information  not  :  sadily 
available  for  devices,  many  were  hesitant  about  using 
electrothermal  models  for  large  scale  testing  of  devices. 

One  study  sponsored  by  the  USAF  attempted  to  apply  a 
much  reduced  version  of  the  electrothermal  models  to  a  wide 
variety  of  devices  where  specific  design  aspects  (such  as 
mask  architecture,  doping  levels  and  mask  areas)  were 
systematically  varied  with  another  aspect  held  as  a  control 
C13).  Results  of  this  large  scale  effort  did  not  show  any 
universal  trend  in  the  power  dissipated  versus  time  to  second 
breakdown  for  any  particular  device  design  aspect.  flany 
individual  results  showed  multiple  peaks  on  graphs  of  number 
failed  versus  power  to  failure  Csee  figures  3  and  4)  and  some 
results  did  not  appear  to  match  the  normal  distribution 
assumed  for  the  probability  of  failure  assessment  of  the 
devices  (many  attempts  to  model  device  failure  sought 
statistical  methods  to  account  for  the  broad  range  of  power 
to  failure  found  experimentally). 

Although  the  variation  in  power  to  failure  had  been  seen 
before,  it  was  believed  to  be  due  to  inter-vendor  differences 
in  device  design  C14II  as  shown  in  figure  5. 
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Figure  5.  Failure  pau/er  threshold  varia¬ 
bility  versus  manufacturers  for  2N*13B9  transis¬ 
tors.  Circles  represent  the  average  value  of 
the  failure  power.  Bars  indicate  the  range  of 
failure  found.  Devices  were  tested  with  the 
emmiter  base  in  reverse  bias  and  the  collector 
open.  From  reference  14. 
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Yet,  the  USAF  sponsored  work  did  not  shoui  any  expected 
correlation  between  device  design  and  the  distribution  of 
failure . 

Critics  of  this  study  believed  that  it  suffered  the  same 
drawbacks  as  previous  experimentation  with  thermal  models; 
namely,  poor  experimental  design  that  did  not  Cor  could  not) 
control  the  primary  parameters  that  effect  failure. 

These  critics,  who  were  investigating  conclusion  two, 
also  pointed  out  that  second  breakdown  and  device  failure  by 
thermal  events  were  not  synonymous.  In  addition,  the  models 
used  were  based  on  thermal  diffusion  and  were  more  applicable 
to  device  failure  due  to  melting,  not  the  initiation  of 
second  breakdown . 

These  critics  had  also  been  able  to  show  experimentally 
that  the  fine  structure  seen  upon  the  initiation  of  second 
breakdown  was  most  certainly  occurring  before  melting  and  was 
most  likely  to  be  due  to  electronic  processes  12,6,7,1511. 

The  most  researched  cause  of  second  breakdown  by  this 
group  involved  avalanche  breakdown  in  conjunction  with 
heating,  defects  or  surface  effects.  These  conjunctive 
causes  were  viewed  as  necessary  to  nucleation  of  the  current 
to  a  current  filament. 

Further  evidence  in  favor  of  this  approach  was  obtained 
by  Sunshine  and  Lampert  C163  that  demonstrated  that  avalanch¬ 
ing  silicon-on-sapphire  diodes  repeatably  showed  multiple 
filament  formation  dependent  on  the  degree  of  electrical 


stress  applied  to  the  device.  By  the  early  19B0s  it  uas 
clear  that  second  hreakdouin  was  the  precursor  of  device 
failure.  However,  no  model  was  able  to  explain  completely 
both  the  characteristic  low  voltage,  high  current  state  and 
the  formation  or  collapse  of  current  to  a  small  filamentCs). 
In  fact,  some  researchers  believed  that  avalanche  breakdown 
and  device  heating  would  not  allow  filament  formation  due  to 
thermal ization  of  the  carrier  distribution.  If  this  was 
found  to  be  true,  no  simple  solution  appeared  to  be  able  to 
explain  the  complex  phenomena  accompanying  second  breakdown. 

Preliminary  Results  of  a  Formal 
Redevelopment  of  Device  Failure 
Phenomena 

After  a  perusal  of  past  research,  the  author  of  this 
work  set  out  to  re-examine  the  problem  of  explaining  device 
behavior  preceding  failure  by  finding  a  physical  explanation 
for  the  current  f i lamentation  seen  under  second  breakdown. 

The  starting  point  for  this  effort  was  basic  carrier 
transport  under  low  and  high  voltages  in  insulators  and 
semiconductors.  These  states  of  electrical  conductivity  in 
materials  were  chosen  because  an  active  semiconductor  device 
is  dominated  by  the  electrical  activity  in  the  depletion 
region,  which  can  exhibit  insulator  or  semiconductor  charac¬ 
teristics  dependent  on  an  external  bias.  Also,  second 
breakdown  in  devices  appeared  exclusively  in  depletion 
regions  as  determined  from  experiments  C17]  and  postmortem 


Explaining  cLirrant  transport  in  a  device  in  this  manner 
will  allow  U3  to  determine  when,  and  if,  a  negative  resist¬ 
ance  regime  has  been  achieved  in  a  sample  or  device. 

To  begin  this  analysis,  we  will  describe  the  basic 
transport  mechanisms  for  space  charge  limited  single  carrier 
injection  in  an  insulator.  This  will  be  followed  by  a  short 
description  of  two  carrier  currents  in  insulators  and 
semiconductors.  Each  case  comes  closer  to  describing  current 
flow  in  the  depletion  region  of  a  device  or  a  sample  of 
material.  These  treatments  follow  from  the  work  of  Hark  and 
Lampert  CIS]  and  were  extensively  used  by  the  author  in  a 
previous  effort  CBO] .  Since  these  derivations  appear 
elsewhere  their  results  will  simply  be  stated  here. 

The  first  case  of  interest  deals  with  single  carrier 
injection  in  insulator-like  materials.  To  approximate  the 
effect  more  closely  to  a  semiconductor  device’s  depletion 
region,  the  insulator  we  discuss  has  connections  that  are 
ohmic  to  electrons  only.  Under  these  conditions,  the 
insulator  is  found  to  transition  through  four  different  types 
of  space  charge  limited  current  flow  with  increasing  applied 
Cdirect  current)  voltage.  At  low  voltages,  current  follows 
Ohm’s  law.  As  the  voltage  increases,  more  electrons  are 
injected  and  begin  to  fill  any  traps  in  the  material  near  the 
cathode.  Yet,  an  Insufficient  number  of  electrons  exist, 
after  trapping,  to  transit  the  material.  Hence,  a  buildup  of 
space  charge  begins  near  the  cathode.  Increasing  the  applied 


voltage  further  ujill  fill  more  traps  and  some  electrons  make 
it  through  the  material  to  the  anode.  Eventually,  the 
voltage  is  raised  enough  so  that  all  traps  in  the  material 
are  filled  and  the  current  follows  a  trap-free  square  law 
characteristic  Csee  figure  6).  Electrons  are  now  able  to 
recombine  primarily  at  the  anode  and  little  trapping  takes 
place  in  the  bulk. 

These  results  are  quite  similar  to  current  flow  in  a 
forward  biased  diode  CSl]  where  current  flow  is  dominated 
primarily  by  recombination  of  minority  carriers  at  the 
depletion  layer  edge  Csee  figure  7).  Whether  or  not  a  diode 
follows  a  trap-free  square  law  is  unclear  at  present. 

If  we  remove  the  constraint  of  ohmic  contacts  for 
electrons  only,  space  charge  build-up  in  the  previously 
modeled  insulator,  is  eliminated.  The  insulator  is  now 
modeled  with  a  cathode  that  injects  electrons  and  an  anode 
that  injects  holes.  The  only  limiting  mechanism  to  current 
flow  in  this  situation  is  recombination.  This  recombination 
occurs  via  two  primary  mechanisms  which  are  dependent  on  the 
applied  potential.  Some  recombination  takes  place  at  traps 
in  the  material  during  low  and  moderate  voltages.  When  the 
applied  voltage  is  large  enough,  carrier  transit  times 
becomes  equal  to  the  lifetime,  through  recombination,  and 
carriers  recombine  primarily  at  the  injecting  contacts. 

The  experimental  results  found  for  the  above  case  of 
ohmic  contacts  for  electrons  and  holes  CEOD  is  somewhat 


Figure  6.  Example  of  the  transitions  for 
single  carrier  space  charge  limited  currents  in 
an  insulator.  Calculated  values  are  from 
reference  19. 
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Figure  7.  Forward  current — voltage 
characteristics  for  alloyed  germanium  diodes. 
Current  is  proportional  to  the  voltage  squared. 
This  indicates,  from  the  previous  discussion, 
that  current  is  dominated  by  single  carrier 
injection  and  recombination  at  the  contacts. 
From  reference  21. 


similar  to  single  carrier  injection  results.  At  low  voltages 
the  current  Follows  a  square  law  very  similar  to  the  trap- 
free  -square  law  relation.  In  essence,  we  have  a  single 
carrier  current  until  a  balance  is  reached  with  the  traps, 
which  preferentially  block  the  motion  of  one  carrier,  until 
all  traps  are  filled.  As  the  voltage  increases,  the  current 
becomes  dominated  by  both  carriers  Cthe  traps  are  now  filled) 
to  produce  a  double  injection  cube  law  for  current  flow. 

We  next  bring  attention  to  a  semiconductor  material 
where,  unlike  an  insulator,  a  number  of  thermally  generated 
free  carriers  are  present.  We  model  this  case  with  an 
electron  injecting  cathode  and  a  hole  injecting  anode  in  a 
semiconductor  that  is  slightly  dominated  by  electrons  or 
holes.  Results  of  calculation  show  the  current  Follows  a 
square  law  which  matches  experimental  results  quite  well  for 
a  low  number  of  carriers  CSE]  Csee  figure  0) .  When  a  large 
number  of  carriers  are  present  (from  injection  or  thermal 
production)  the  current  follows  a  cubed  law  CE33  Csee  Figure 
9) .  An  examination  of  the  carrier  concentration  in  both 
cases  of  injection  indicates  a  build-up  of  carriers  at  one  of 
the  material  contacts  Cdependent  on  the  dominant  carrier 
sign)  is  occurring.  This  effect  creates  an  effective  anode 
or  cathode  which  makes  the  semiconductor  seem  electrically 
shorter.  In  a  sense,  a  portion  of  the  sample  experiences  a 
space  charge  limited  current  Flow,  and  the  rest  of  the  sample 
experiences  a  trap  free  current  characteristic. 
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Figure  B.  Current  density  versus  voltage 
characteristics  for  double  injection  into 
lightly  doped  p-type  silicon.  Note  the  square 
iauj  regime.  From  reference  BE. 


WiCSWv’WW 


J 


VOLTAGE,  Y 


Figure  3.  Current  density  versus  reverse 
voltage  characteristics  for  PIN  diodes  at  BOOK, 
Note  the  expected  cube  law  regime  holds  for  a 
device  also.  From  reference  23. 
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We  nou)  deal  with  a  semiconductor  with  contacts  as 
discussed  before,  but,  we  add  a  set  of  recombination  centers 
deep  in  the  forbidden  zone.  This  is  usually  accomplished  by 
placing  additional  dopants  in  the  material  during  the 
crystallization  process  to  provide  an  additional  set  of  traps 
for  recombination.  These  traps  are  usually  formed  in  the 
forbidden  zone  a  fraction  of  an  electron  volt  from  the  energy 
bands.  At  room  temperature,  these  levels  are  empty. 

Analysis  shows  the  presence  of  these  shallow  recombination 
centers,  presents  an  additional  potential  that  the  minority 
carriers  must  uvercome  with  the  aid  of  an  external  voltage. 
Until  the  externally  applied  voltage  is  greater  than  this 
recombination  potential  current  follows  a  trap-free  space 
charge  limited  square  law.  Once  the  recombination  barrier  is 
surpassed  by  the  applied  voltage,  the  sample  potential  drops 
precipitously  and  follows  a  square  law  as  arrived  at  earlier 
for  a  semiconductor  C843  Csee  figure  10). 

When  the  recombination  centers,  added  to  the  sample  as 
doping  during  crystallization  process,  form  traps  near  the 
center  of  the  forbidden  zone  they  exist  in  a  partially  filled 
state  at  room  temperature.  Analysis  shows  that  both  carriers 
now  see  a  potential  barrier  that  must  be  overcome  by  the 
external  voltage  before  current  can  flow  II853  Csee  figure 
11)  . 

These  last  two  situations  with  semiconductors  show 
current  and  voltage  characteristics  which  exhibit  signs  of 


Figure  10.  General  current  density  versus 
voltage  for  the  case  of  double  injection  in  a 
semiconductor  with  recombination  centers. 

Ratios  of  the  carrier  capture  cross-sections  are 
listed  for  electrons  and  holes.  Calculated  by 
Lampert  and  Hark  using  general  parameters.  From 
reference  13. 
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Figure  11.  General  current  density  versus 
voltage  for  double  injection  in  a  semiconductor 
uiith  partially  filled  recombination  centers. 

The  transition  to  negative  resistance  occurs  at 
a  loujer  current  density  than  seen  in  figure  10. 
Calculated  by  Lampert  and  Mark  in  reference  13. 
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what  is  called  negative  resistance.  Negative  resistance,  in 
the  above  cases,  exhibits  a  transition  from  a  high  voltage, 
low  current,  to  a  low  voltage,  high  current  state  quite  like 
the  characteristics  seen  during  second  breakdown.  Since  the 
previous  analysis  for  single  and  double  carrier  currents  is 
applicable  to  the  depletion  region  of  devices,  it  can  be 
inferred  that  second  breakdown  in  devices  is  the  transition 
from  single  carrier  to  double  carrier  injection,  creating  a 
negative  resistance  characteristic. 

Ridley  CS6I1  demonstrated  that  current  controlled 
negative  resistance  CCCNR)  creates  a  regime  where  current 
flow  is  unstable  against  the  formation  of  a  current  filament. 
This  type  of  negative  resistance  is  called  CCNR  after  the 
multi-valued  current  that  can  be  exhibited  for  a  single 
voltage.  In  a  similar  fashion,  Ridley  LEBD  demonstrated  that 
a  under  a  voltage  controlled  negative  resistance  CVCNR), 
current  is  unstable  against  formation  of  current  domains. 

VCNR  is  named  after  the  multi-valued  voltages  that  occur  for 
a  particular  current  level.  Both  aspects  of  negative 
resistance  have  been  verified  experimentally  to  lead  to 
filament  formation  under  CCNR  CB7,E8D  and  domain  formation 
under  VCNR  CE3,30]. 

If  second  breakdown  is  the  initiation  of  CCNR,  then 
second  breakdown  could  be  analyzed  through  the  device’s 
material  parameters.  Yet,  Ridley’s  proof  of  filament 
formation  is  macroscopic,  from  a  thermodynamic  point  of  view. 
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and  a  further  explanation  would  be  required  in  order  to 
characterize  the  filament  in  terms  of  a  model. 

Only  one  physical  mechanism  appeared  capable  of  produc¬ 
ing  the  filamentary  current  channel  seen  in  experiments: 
self-induced  magnetic  pinch  of  the  current  column.  Initial¬ 
ly,  this  premise  was  considered  quite  controversial.  Yet, 
research  in  plasma  physics  applied  to  semiconductors  has 
shown  modeling  current  flow  as  a  solid  state  plasma  is  quite 
realistic . 

Parallels  Between  Classical  Plasmas 
and  Solid  State  Plasmas 

To  demonstrate  the  utility  of  modeling  current  flow  in  a 
semiconductor  as  a  solid  state  plasma,  a  short  history  of  the 
development  of  plasma  physics  and  its  analogous  counterparts 
in  the  development  of  solid  state  plasma  physics  will  be 
examined . 

The  results  of  this  section  will  lay  the  foundation  for 
the  approach  in  later  chapters  to  model  the  solid  state 
plasma  in  terms  of  well  known  classical  plasma  parameters. 

The  term  plasma  was  first  applied  by  Langmuir  C31I1,  in 
19S0,  to  the  ionized  gas  in  an  electrical  discharge.  A 
plasma  was  dominated  by  coulomb  forces  between  individual 
parts  of  the  plasma  that  gave  it  a  cohesiveness  that  has  been 
termed,  collective  behavior.  Following  these  discoveries 
Bennett  C3S1,  in  1934,  predicted  the  self-induced  compression 
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of  an  electron  beam  could  be  produced  by  its  self-magnetic 
field.  Tonks  C331  called  this  behavior  the  pinch  effect. 

Later  in  1335,  Eckersiey  C341  determined  the  dispersion 
relationship  for  electromagnetic  waves  that  propagate  through 
the  ionosphere  along  the  Earth’s  magnetic  field  lines.  The 
group  velocity  of  these  waves  increases  with  frequency  such 
that  an  observer  detecting  such  waves  finds  a  tone  that  falls 
with  time.  Barkhausen  C351  had  detected  these  waves  as  early 
as  Ulorld  War  I  and  coined  the  term  by  which  they  are  pre¬ 
sently  called:  whistlers. 

In  134S  Alfvdn  C3B1  predicted  a  new  type  of  electromag¬ 
netic  wave  propagating  through  a  highly  conductive  plasma  in 
a  strong  magnetic  field.  AlFvfen  determined  this  from  a  study 
of  the  phenomena  that  occur  near  sunspots.  Their  subsequent 
discovery  in  the  laboratory  has  lead  to  the  appropriate  name: 
Alfvdn  waves. 

Plasmas  in  solids  were  first  postulated  by  the  studies 
of  Rutheman  C371  and  Lang  C381  in  1940.  In  their  experi¬ 
ments,  the  passage  of  an  electron  beam  through  a  metal  foil 
exhibited  energy  losses  that  were  discrete  in  nature.  Later, 
in  195E,  this  effect  was  expiained  by  Bohm  and  Pines  11331  as 
the  excitation  of  electron  plasma  oscillations  in  the  metal 
in  the  form  of  plasmons  Ca  quantized  oscillation  of  the 
electron  plasma) . 

In  the  late  1350 ’s  and  early  1960 ’s,  a  number  of  results 
were  reported  concerning  anomalous  electrical  behavior  in 
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germanium,  indium  antimonide,  gallium  arsenide  and  bismuth. 
Oscillations  were  reported  in  indium  antimonide  and  gallium 
arsenide  under  conditions  of  impact  ionization  and  in  the 
presence  of  a  longitudinal  magnetic  field  that  was  comparable 
in  magnitude  to  the  intrinsic  magnetic  field  at  the  samples’ 
surface  C40,411.  Studies  of  germanium  at  lou  temperatures 
and  strong  magnetic  fields  shouted  lout  frequency  current 
oscillations  CiSl .  Absorption  of  electromagnetic  waves  in 
the  centimeter  range  in  bismuth,  under  strong  magnetic 
fields,  exhibited  a  linear  dependence  of  the  absorption 
coefficient  with  magnetic  field  strength  114311. 

Researchers  found  that  these  phenomena  could  not  be 
explained  in  terms  of  simple  current  transport  variables  and 
a  new  approach  certainly  seemed  warranted. 

The  results  found  for  germanium  under  high  magnetic 
fields  were  soon  explained  by  Glicksman  C44I]  on  the  basis  of 
treating  current  flow  as  a  solid  state  plasma  using  the 
theory  of  the  helicon  instability  Csimiiar  to  whistlers  in 
the  ionosphere) .  The  results  in  bismuth  were  determined  to 
be  due  to  propagation  of  an  AifvPn  wave  whose  phase  velocity 
is  proportional  to  the  magnetic  field  strength  which  governs 
the  absorption  coefficient  C45I1  . 

Finally,  the  results  in  indium  antimonide  and  gallium 
arsenide  were  explained  by  Glicksman  and  Steele  C45,4BI1  with 
the  use  of  the  pinch  effect.  In  their  work,  an  electron-hole 
plasma  is  strongly  compressed  by  the  sample’s  current 


interaction  with  its  self-magnetic  field.  As  the  density 
rises  in  the  column,  the  resistance  increases  as  enhanced 


electron-hole  scattering  and  bulk  recombination  occur.  The 
study  also  demonstrated  that  the  criterion  for  compression 
developed  by  Bennett  C351  was  well  satisfied  experimentally 
C47] .  Later,  Glicksman  and  Ando  C483  demonstrated  that  under 
a  strong  enough  longitudinal  magnetic  field,  the  pinch 
transformed  into  a  helicon  instability.  In  1363,  Gunn  II433 
determined  that  micrawave  radiation  was  generated  by  gallium 
arsenide  and  indium  phosphide  under  strong  electric  fields. 
The  effect,  which  was  named  after  Gunn,  was  found  to  be  due 
to  VCNR  and  was  predicted  earlier  by  Ridley  CE63  .  The  Gunn 
effect  had  no  counterpart  in  gaseous  plasma  physics  experi¬ 
ments.  Yet,  a  rigorous  derivation  of  the  Gunn  effect  was 
carried  out  using  the  formalism  of  collisionless  shock  waves 
developed  by  Sagdeev  [503.  Later,  an  effect  analogous  to  the 
Gunn  effect  was  found  in  gaseous  plasmas  in  the  study  of  ion 
acoustic  shock  waves  C513.  Here  was  a  case  where  develop¬ 
ments  in  solid  state  plasmas  were  a  precursor  to  discovery  of 
the  same  effects  in  gaseous  plasma  experiments. 

At  this  point,  it  should  be  quite  clear  to  the  reader 
that  although  there  are  differences  between  a  solid  state 
plasma  and  a  gaseous  plasma,  the  approach  to  problems  and 
much  of  the  phenomena  they  share  Csuch  as  instabilities)  are 
grounds  of  commonality,  riodeling  current  flow  in  a  semicon¬ 
ductor  as  a  solid  state  plasma  is  well  qualified. 


If  the  reader  accepts  the  previous  argument  by  analogy 
of  the  correspondence  between  gaseous  plasmas  and  solid  state 
plasmas,  then  second  breakdown  in  semiconductor  devices  could 
be  explained  by  modeling  current  flow  as  a  solid  state  plasma 
under  CCNR  conditions.  Fi lamentation  could  then  be  explained 
as  a  result  of  the  plasma  pinch. 

The  balance  of  this  work  will  be  devoted  to  determining 
if  pinching  of  the  solid  state  plasma  does  occur  under  CCNR 
conditions,  whether  or  not  the  pinch  is  stable,  and  what 
additional  effects  could  ensue  if  the  pinch  forms. 

Chapter  II  will  start  with  a  detailed  examination  of 
plasmas  to  determine  a  set  of  parameters  to  characterize  a 
solid  state  plasma  in  a  way  that  is  analogous  to  known 
parameters  for  a  classical  plasma. 

Chapter  III  will  examine  the  conditions,  and  limita¬ 
tions,  CCNR  places  on  the  environment  that  contains  the  solid 
state  plasma. 

Chapter  IV  will  model  a  solid  state  plasma  under  CCNR 
conditions  and  determine  if  pinching  occurs,  and  whether  or 
not  it  is  stable.  This  analysis  will  develop  a  set  of 
criteria  to  determine  the  initiating  conditions  for  the 
pinch . 

Chapter  V  will  look  at  experimental  measurements  of  a 
transistor  that  was  driven  electrically  to  the  initiation  of 
second  breakdown,  then  recovered  safely.  Analysis  of  the 
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device  physics  uill  be  performed  to  determine  if  pinching 


conditions  were  achieved  at  the  initiation  of  second  break¬ 


down  . 
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CHAPTER  II 


( 


GENERAL  CHARACTERIZATION  OF  THE  SOLID  STATE 


PLASMA:  PLASMA  CRITERIA 


It  was  praposed  in  chapter  I  that  the  electron-hole 
cloud  in  a  semiconductor  could  be  treated  as  a  solid  state 
plasma.  Parallels  mere  drawn  between  developments  in  clas¬ 
sical  plasma  physics  and  experimental  work  that  demonstrated 
similar  plasma  characteristics  in  solids. 

This  chapter  will  elaborate  on  the  connection  between 
classical  plasmas  and  solid  state  plasmas  by  utilizing  a  set 
of  parameters  to  characterize  all  plasmas. 

The  well  known  classical  plasma  parameters  Ce.g.,  Debye 
length  and  Debye  sphere!  will  be  Found  to  represent  a  subset 
of  the  most  general  parameters  a  plasma  can  possess. 

Application  of  these  new  parameters  to  a  solid  will  show 
the  solid  state  plasma,  although  quantum  mechanical  in 
origin,  can  be  treated  as  a  weakly-coupled  classical  plasma 
at  room  temperature,  in  intrinsic  germanium,  silicon  and 
gallium  arsenide. 


The  Classical  Plasma 

A  plasma  is  generally  defined  as  an  ionized  gas  of 


electronically  quasi-neutral  particles  that  exhibit  a 


collective  behavior  in  reaction  to  the  surrounding  environ¬ 
ment  C52].  The  term  quasi-neutral  refers  to  a  neutrality 
that  occurs  due  to  the  kinetic  action  of  the  particles 
forming  an  electrostatic  shield  batueen  unlike  charges.  As  a 
result  of  this  shielding,  the  potential  field  from  a  test 
charge  placed  in  the  plasma  falls  off  faster  than  in  a 
vacuum.  This  action  is  also  a  direct  consequence  of  the 
kinetic  energy  of  the  particles  being  greater  than  the 
potential  energy  between  particles.  By  solving  Poisson’s 
equation  for  a  plasma,  or  looking  at  the  ratio  of  kinetic  and 
potential  energies  of  the  particles,  a  simple  relationship 
evolves  specifying  a  distance  Csimilar  to  a  skin  depth)  over 
which  external  potentials  vanish.  This  action  is  known  as 
Debye  shielding  and  forms  a  critical  concept  necessary  to  the 
definition  of  a  plasma.  The  term  collective  behavior 
describes  how  a  plasma  reacts  as  a  single  system  to  the 
surrounding  environment  in  order  to  maintain  its  shielding 
characteristics . 

In  order  to  differentiate  between  a  plasma  and  a 
collection  of  free  particles,  a  set  of  criterion  have  been 
derived,  based  on  classical  grounds,  to  determine  the  degree 
of  collective  behavior  in  the  presence  of  electromagnetic 
f ields , 

The  simplest  criterion  deals  with  the  Debye  length.  If 
the  entire  system  of  charges  is  larger  than  the  Debye  length 
then  Debye  screening  of  external  electromagnetic  fields  is 
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meaningful.  In  CGS  units  the  Debye  length  takas  the  form 


XjD  = 


KeTe 

4Tme^ 


<C  the  systems ’s  dimensions  , 


where 


Kb  •  Boltzmann’s  constant. 

T  -  The  carrier  temperature, 
n  -  The  carrier  density, 
e  “  The  electronic  charge, 
e  “  The  material’s  dielectric  constant. 

A  second  criterion  determines  if  there  are  enough  charge 
carriers  within  a  sphere  of  radius  equal  to  the  Debye  length 
to  accomplish  shielding  of  external  fields.  This  criterion 
can  be  derived  from  the  carrier  density  and  the  definition  of 
a  sphere  to  arrive  at 


Nb  —  —mrX^  ^  1  electron  . 


The  final  criterion  deals  with  a  plasma’s  ability  to 
possess  collective  behavior  that  is  not  swamped  out  by  its 
own  random  collisions.  To  fulfill  this  criterion,  the 
product  of  the  frequency  of  oscillation  in  the  plasma  C  w  ) 
and  the  mean  time  between  collisions  C  r  )  must  be  greater 
than  unity .  This  relationship  can  be  written  as 
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These  criteria  C  Ad  ,  Nd  and  tor  ),  also  detail 
general  characteristics  of  a  plasma  that  are  useful  in 
dealing  with  the  scale  and  types  of  flux  and  forces  at  work 
in  a  plasma.  It  is  easy  to  graph  these  relations  to  see  the 
general  scope  that  plasmas  cover  using  these  classical 
criteria.  Figure  IS  clearly  shous  known  plasmas  exist  for  a 
broad  range  of  temperatures  and  carrier  densities.  It  is 
also  clear  that  the  dominant  forces  driving  any  particular 
plasma  will  differ  depending  on  the  plasma’s  characteristics. 
For  example,  at  low  temperatures  coulomb  forces  are  important 
while  at  high  temperatures  thermal  forces  dominate.  At  high 
densities  quantum  mechanical  forces  cannot  be  ignored,  while 
at  low  densities  carriers  act  as  a  free  classical  gas. 

Since  the  plasma  criteria  discussed  above  are  for  a 
classical  plasma  Ci.e.,  a  gas  that  follows  naxwell-Boltzmann 
statistics)  we  cannot  apply  them  directly  to  plasmas  where 
quantum  mechanical  effects  become  dominant,  as  they  would  in 
a  solid  state  plasma.  Thus,  we  will  have  to  derive  new 
criteria  for  the  solid  state  plasma. 

The  Solid  State  Plasma  Criteria 
The  primary  distinction  between  a  classical  gas  and  a 
quantum  mechanical  gas  is  the  distribution  of  carriers  with 


*• 


FigLire  IE.  Plots  of  classical  plasma 
parameters.  Part  a  is  a  plot  of  the  Debye 
length  in  meters  versus  the  concentration. 
Part  b  is  a  plot  of  the  concentration  within 
the  Debye  sphere  versus  concentration. 


the  system’s  possible  energy  states.  A  classical  gas 
reflects  a  riaxujel l-Bol tzmann  distribution  while  a  quantum 
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mechanical  gas  reflects  a  Farmi-Dirac  distribution  Csince  we 
are  dealing  with  electrons  as  the  prime  carriers  within  a 
solid)  . 

This  distinction  is  important,  since  the  classical  gas 
was  the  basis  for  the  previously  reviewed  plasma  criteria. 
Simple  re-application  of  those  criteria  is  meaningful,  if  the 
carriers  present  in  a  material  can  be  treated  using  Maxwell- 
Boltzmann  statistics.  If  conditions  do  not  exist  to  favor  a 
riaxwel  1-Boltzmann  distributjon  the  gas  may  be  affected  by 
quantum  mechanical  effects,  and  a  quantum  mechanical  distri¬ 
bution  would  be  required. 

To  accommodate  this  possibility  it  is  necessary  to  be 
able  to  distinguish  between  plasmas  utilizing  a  classical  or 
a  quantum  mechanical  approach  r53,541.  The  simplest  dividing 
post  is  a  comparison  of  the  chemical  potential  or  Fermi 
energy  and  the  thermal  energy.  Classically,  we  have  the 
kinetic  energy 


The  De  Brogile  relation  allows  us  to  rewrite 


equation  4  to  arrive  at 
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Ef  =  -^~- 


Where  Ef  is  the  Fermi  energy.  The  Fermi  energy  generally 
represents  the  energy  of  carriers  from  a  quantum  mechanical 
perspective.  Comparison  of  the  Fermi  energy  and  a  classical 
kinetic  energy  ujould  determine  which  perspective  dominates 
carrier  activity .  A  simple  calculation  shows  that  the  Fermi 
energy  for  copper  is  around  6.9  eV  compared  to  the  thermal 
energy  at  room  temperature  of  about  O.OEl  eV .  Quantum 
mechanical  effects  certainly  cannot  be  ignored  since  the 
energy  from  a  quantum  mechanical  perspective  is  over  300 
times  greater  than  the  classical  kinetic  energy.  If  we 
express  the  Fermi  energy,  via  the  De  Brogile  relationship, 
as  a  temperature,  we  arrive  at  the  Familiar  equation  for  the 
De  Brogile  wavelength 


2Trm’'  K  qT 


Where 


771*  =  The  effective  mass 


T  =  The  carrier  temperature 


The  effective  mass  reflects  the  quantum  mechanical 
effect  of  the  carrier  moving  through  the  lattice  potentials, 
which  alter  the  carrier’s  response  to  an  external  electromag¬ 
netic  Field. 
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If  the  average  distance  between  carriers  is  on  the  order 
of,  or  smaller  than  Xqq  ,  the  gas  is  quantum  mechanical  in 
behavior.  In  explicit  terns  we  can  define  this  as 

C7) 

If  A  is  greater  than  or  equal  to  unity,  the  gas  is 
quantum  mechanical  and  if  A  is  much  smaller  than  unity,  it 
behaves  like  a  classical  gas.  Lise  of  the  De  Brogile  wavelen¬ 
gth  gives  a  well  known  rule  of  thumb  measure  of  the  quantum 
mechanical  effects  that  could  occur  in  various  situations. 
Application  of  the  relationship  shows  quantum  mechanical 
effects  dominate  at  high  dcjnsities  and  low  temperatures.  At 
low  densities  and  high  temperatures,  quantum  mechanical 
effects  are  minimized.  This  criterion  allows  a  distinction 
to  be  made  between  a  classical  plasma  and  a  quantum  mech¬ 
anical  plasma. 

Our  Oebye  sphere  Cequatinn  requirement  could  be 
rewritten  as 

Nd  --=  ^’rnApB  "=  CBl 

where  F  must  be  less  than  unity  to  satisfy  the  criterion 
for  a  classical  plasma,  Ini.erent  in  this  derivation  is  the 
understanding  that  the  pater.rial  energy  between  two  carriers 
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is  small  compared  to  their  kinetic  energy.  This  can  be  seen 
by  taking  the  ratio  of  the  two  energies  as  folloujs  Cin  CBS 
units ) 


potential  energy  2e^n3 
kinetic  energy  ZeK^T' 


potential  energy 
kinetic  energy 


CIO) 


Where  C  is  some  constant  and  e  is  the  material’s  dielectric 
constant . 

The  new  plasma  parameter  F  reflects  the  fact  that  the 
potential  energy  between  particles  should  be  much  smaller 
than  the  kinetic  energy  of  a  particle  to  satisfy  our  previous 
criteria  for  classical  plasma  behavior  C  F  <<  1  implies 
Nd  1  ) .  A  plasma  with  a  F  much  smaller  than  unity  is 
called  a  weakly-coupled  plasma  C55D . 

The  case  of  F  greater  than  or  equal  to  unity  corres¬ 
ponds  to  what  is  called  a  strongly-coupled  plasma.  This  case 
does  not  have  a  classical  analog  since  this  type  of  plasma  is 
controlled  by  the  coulombic  potential  between  individual 
particles.  Inspection  of  equations  7  and  0  for  a  strongly- 
coupled  plasma  indicates  a  very  high  carrier  density  must  be 
present  11563  and  quantum  mecnanical  forces  will  dominate. 


Very  high  carrier  densities  in  a  solid  would  see  quantum 
mechanical  affects  from  the  lattice  potentials  and  a  spread 
in  carrier  momenturns  due  tc  the  Heisenburg  Uncertainty 
Principle.  We  can  quantify  these  affects  on  the  kinetic 
energy  as 


(Ap)^  ^  h~  1  h? 
2m  2m*  2m* 


Equation  11  can  be  restated  in  terms  of  the  Fermi  energy  as 


To  satisfy  our  definition  of  a  plasma  the  potential  energy 
between  particles  must  be  lass  than  the  kinetic  energy . 
Taking  the  ratio  of  potentiaJ  to  kinetic  energies  yields 


4-  (3 


c  =  _  _ 

^/'e(37r2)2/3^l/3• 


The  ratio  Cp  is  equivalent  to  the  inverse  of  the  Debye 
sphere  in  a  classical  plasma,  but,  since  the  plasma  is 
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quantum  mechanical  in  origin,  it  is  called  a  degenerate 
plasma.  If  Cp  is  less  than  or  equal  to  unity,  the  plasma  is 
a  weakly-coupled  degenerate  plasma  (analogous  to  No  >  ^  in  a 
classical  plasmaD  .  If  Cp  is  greater  than  unity  then  the 
plasma  is  a  strongly-coupled  degenerate  plasma  (which  has  no 
classical  counterpart) .  Cp  can  also  be  arrived  at  by 
comparing  the  radius  of  the  effective  Bohr  orbit  (for  n  ”  1) 
to  the  interparticle  spacing  as. 
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Here,  r,  is  related  to  C-p  by 


J 


Liboff  C57D  calls  r,.  the  quantum  compression  parameter. 
Values  of  r,.  less  than  or  equal  to  unity  correspond  to  a 
uieakly-coupled  degenerate  plasma  and  values  greater  than 
unity  corresponds  to  a  strongly-coupled  degenerate  plasma. 
When  r,  is  greater  than  unity  the  plasma  approaches  a 
quantized  state  Cplasmons  are  formed).  Similar  statements 
can  be  made  for  Cp  • 

We  can  now  examine  classical  and  degenerate  plasmas  in 
terms  of  the  new  plasma  parameters. 

As  shown  in  figure  13,  classical  plasmas  correspond  to  A 
and  r  much  smaller  than  unity,  that  is,  they  are  weakly- 
coupled  classical  plasmas.  Degenerate,  or  quantum  mechanical 
plasmas,  are  characterized  by  A  equal  to,  or  greater  than 
unity.  Application  of  r.  Cor  Cp  )  indicates  the  degree  of 
quantization  the  degenerate  plasma  possesses.  A  value  of  r, 
greater  than  unity  indicates  the  plasma  is  highly  quantized 
Cplasmons  form,  equivalent  to  a  strongly-coupled  plasma),  and 
values  of  r,  less  than  unity  indicate  the  plasma  is  approach¬ 
ing  the  classical  characteristics  of  a  weakly-coupled  plasma. 

Redrawing  the  piots  of  figure  IS  in  terms  of  the  quantum 
degeneracy  parameter,  A  , the  plasma  coupling  parameter,  r  i 
and  the  quantum  compression  parameter,  r,  indicates  how  the 
carrier  distributions  for  a  system  of  particles  dictates  the 
type  of  plasma  behavior  obtained. 

Figure  14  shows  a  plot  of  concentration  versus  the 
quantum  degeneracy  parameter.  This  plot  shows  how  specific 
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Figure  13.  Relationships  For  the  neiu 
plasma  parameters.  Thiese  relationships  are 
more  general  than  chose  utilised  for  a  clas 
sical  plasma  alone. 
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details,  such  as  the  ccncentration  and  temperature,  result  in 
either  a  classical  or  degenerate  plasma  .  The  quantum 
compression  parameter  shows  how  solid  state  plasmas  vary  from 
a  highly  quantized  plasma  in  a  metal,  to  a  classical  plasma 
in  some  semiconductors. 

Figure  15  shows  a  plot  of  the  carrier  concentration 
versus  the  plasma  coupling  parameter.  This  plot  indicates 
the  close  analogy  between  and  F  as  they  are  applied  to  a 

degenerate  or  classical  plasma.  For  a  solid  state  plasma 
with  low  carrier  concentrations  we  can  utilize  classical 
results,  since  this  area  corresponds  to  a  weakly-coupled 
plasma . 

These  graphs  indicate  that  in  order  to  call  a  collection 
of  particles  a  plasma  they  must  exhibit  collective  behavior 
in  some  surrounding  environment.  For  example,  an  interstel¬ 
lar  gas,  such  as  the  solar  wind,  may  have  a  concentration  of 
5  electrons  per  cubic  centimeter.  With  an  electron  tempera¬ 
ture  around  50  eV  and  an  ion  temperature  near  10  eV  the  Debye 
length  for  the  solar  wind  is  around  24  meters  and  the  Debye 
sphere  criteria  is  satisfied.  Yet,  the  solar  wind  does  not 
qualify  as  a  classical  plasma  until  we  talk  of  scales  much 
larger  than  24  meters.  At  these  scales,  the  kinetic  energy 
of  the  particles  are  much  larger  than  the  potential  energy  of 
interaction,  which  satisfies  our  definition  of  a  plasma. 

Similar  reasoning  applies  to  a  solid  state  plasma. 
Calculations  of  the  Debye  length  for  electrons  in  a 


Figure  15.  Plot  of  the  concentration 
versus  the  plasma  parameter  for  classical 
coupling.  Some  areas  where  known  plasmas  may 
exist  are  indicated.  The  key  is  as  follows.- 
solar  corona,  b  =  gaseous  discharge,  c  = 
interstellar  gas,  d  =  semiconductors,  e  = 
thermonuclear  plasma,  f  =  metals,  g  =  Jovian 
interior,  h  =  red  giant  star  interior. 
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semiconductor  yield  values  that  are  small  compared  to  a 
laboratory  sample  scale.  But,  such  samples  barely  satisfy 
the  Debye  sphere  criterion.  Large  carrier  concentrations 
force  us  to  account  for  the  quantum  mechanical  source  of  the 
plasma.  That  source  is  the  lattice.  At  values  of  the 
quantum  compression  parameter  near  unity  the  effect  of  the 
ion  potentials  of  the  lattice  are  nil  on  the  plasma,  and  we 
have  a  uieakly-coupled  plasma.  At  larger  values  of  the 
compression  parameter,  the  plasma  is  strongly-coupled,  or, 
the  presence  of  the  lattice  determines  the  plasma’s  charac¬ 
ter.  In  other  words,  the  solid  state  plasma  is  a  plasma 
because  of  the  lattice.  Remove  the  lattice,  or  speak  of 
scales  very  large  compared  to  the  solid  state  plasma,  and 
plasma  behavior  disappears. 

With  these  ideas  in  mind,  similar  graphs  can  be  produced 
for  the  semiconductors  of  interest  in  this  study.  Those 
semiconductors  are;  germanium,  silicon  and  gallium  arsenide. 
Results  are  shown  on  figures  16,  17  and  18  respectively. 

These  semiconductors  are  used  often  in  semiconductor  devices 
and  graphs  of  the  quantum  degeneracy  parameter  versus 
concentration  shows  what  plasma  characteristics  we  could 
expect . 

These  graphs  demonstrate  that  most  of  the  plasma  regime 
for  these  semiconductors  is  classical  in  nature.  This  result 
is  quite  consistent  with  the  Fermi-Dirac  distribution  these 
semiconductors  follow. 


Figure  17.  Plat  of  the  concentration 
versus  the  quantum  degeneracy  parameter  for 
silicon , 
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Figure  18.  Plot  of  the  concentration 
versus  the  quantum  degeneracy  parameter  for 
gallium  arsenide. 
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Due  to  the  structure  of  this  distribution,  only  carriers  in 
the  tail  of  the  distribution  contribute  to  the  semiconduc¬ 
tor’s  electrical  properties  for  temperatures  greater  than 
zero.  Yet,  the  tail  of  the  Fermi-Dirac  distribution  falls 
off  approximately  like  a  (laxuiell-Bol tzmann  distribution. 
Therefore,  the  conclusions,  showing  classical  plasma  behav¬ 
ior,  drawn  are  is  not  in  error. 

It  is  also  interesting  to  observe  the  effect,  in  figures 
16,  17  and  16,  of  the  effective  mass  and  the  dielectric 
constant  on  solid  state  plasma  behavior.  In  most  semicon¬ 
ductors  the  effective  mass  is  usually  less  than  unity  and  the 
dielectric  constant  is  greater  than  unity .  The  combined 
effect  of  these  variables  has  a  general  tendency  to  move  the 
solid  state  plasma  in  semiconductors  toward  greater  degener¬ 
acy  Cor  metallic  behavior)  at  low  temperatures  or  large 
concentrations . 

By  analogy  with  a  classical  plasma,  we  would  also  expect 
a  classically  treated  solid  state  plasma  to  be  most  sensitive 
to  large  scale  perturbations  Con  the  scale  of  the  plasma  in 
the  sample)  .  Highly  degenerate  solid  state  plasmas  would  be 
more  sensitive  to  variations  in  the  distribution  of  carriers. 

Instabilities  at  the  sample  size  could  be  modeled  using 
macroscopic  approaches,  like  magnetohydrodynamics.  Instabil¬ 
ities  at  the  carrier  distribution  level  would  require 
modeling  using  kinetic  theory  via  the  Boltzmann  equation. 


This  analogy  betuieen  classical  and  solid  state  plasmas 


will  allow  us  to  model  solid  state  plasmas  in  close  cor¬ 
respondence  to  weakly  coupled  gaseous  plasmas.  This  will 
certainly  be  true  for  hydrodynamic  instabilities.  Kinetic 
modeling  may  be  similar,  but,  care  would  have  to  be  exercised 
since  the  recombination-generation  mechanisms  present,  in 
either  case,  are  not  analogous. 


chafte;r  III 


ANALYSIS  OF  SPECIFIC  NEGATIVE  RESISTANCE 
IN  SOLIDS  PROPOSED  BY  RIDLEY 

B.  K.  Ridley’s  paper  126]  on  specific  negative  resis¬ 
tance  is  central  to  understanding  the  conditions  surrounding 
CCNR .  Ridley  utilized  an  approach  based  on  arguments  from 
irreversible  thermodynamics,  applied  to  a  situation  where  a 
general  CCNR  characteristic  occurs.  Knowledge  of  the  general 
tendencies  that  occur  in  CCNR  could  then  be  applied  to  second 
breakdown . 

This  chapter  will  examine  the  basic  equations  of  motion 
in  irreversible  processes  to  allow  a  quantification  of  the 
CCNR  state  in  terms  of  specific  variables  such  as  the 
current,  electric  and  magnetic  fields,  and  thermal  farces. 
This  information  is  necessary  in  order  to  obtain  boundary 
conditions  to  model  current  flow  as  a  solid  state  plasma. 

Ridley’s  arguments  will  also  be  examined  in  the  light  of 
criticisms  put  forth  by  Volkov  and  Kogan  C59] ,  and  a  later 
paper  by  Takeyama  and  Kitahara  CBO],  which  casts  doubts  on 
the  validity  of  Ridley’s  approach  compared  to  the  previous 
work  of  De  Groat  and  Prigogine. 

De  Groat’s  derivation  will  he  shown  first  and  compared 
to  Ridley’s  derivation. 


The  state  variables  of  a  continuous  system  can  depend  on 


time  and  space.  For  a  semiconductor,  this  is  certainly  true 
if  temperatures  Cof  the  lattice  and  carriers)  are  alloued  to 
vary  due  to  external  driving  forces,  such  as,  the  electric 
field  and  current.  These  various  driving  forces  and  fluxes 
are  fundamental  variables  relating  conservation  of  mass, 
energy,  balance  of  forces  and  Gibb’s  equation.  The  usual 
forms  for  these  relationships  are 


dt 


V  • {pv)  =  0, 


C19) 


—  (-pV^-\-pU)  -  F  -V  -  V  ■  [Pv  +  Jg) 
DV2 


CEO) 


dv 


P— 


^  -VPPF, 


CEl ) 


1 
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aihere  the  variables  are  identified  as 
P  =  mass  density, 

V  =  velocity, 

=  convective  derivative, 

=  heat  floui, 
p  =  external  forces, 

U  =  internal  energy  Cper  unit  volume), 

V  =  volume, 

T  temperature, 

S  =  entropyCper  unit  volume), 
p  =  pressure, 

/ifc  =  chemical  potential  for  species  k. 

The  general  procedure  relating  the  various  driving 
forces  and  their  resultant  fluxes  has  been  derived  by  De 
GrootC67D  and  is  repeated  here.  De  Groot’s  goal  is  to  arrive 
at  a  continuity  equation  for  the  entropy  flows.  If  a 
stationary  state  Ca  non — equilibrium  state  of  an  open  system) 
exists,  the  principle  of  least  entropy  production  can  be 
applied  to  the  entropy  flows  to  determine  the  general 
direction  the  state  variables  will  take. 

The  law  of  conservation  of  mass  is  written  as 


dpk 

dt 


+  V  •  {pkvu)  = 


1  dMk 

v~ir 


C23) 


rjn 


I 


5^ 


N 


luhere 


J(.  =  chemical  reaction  rate  in  masc  per  unit  '.'nlume  and 
time  . 

Pfc  =  mass  density  of  the  kth  component, 

=  velocity  of  the  kth  component, 

Mk  =  mass  of  the  kth  component, 

=  stoichiometric  number  of  components. 


d 

I 


■*  «><•  i’l 


k**  1*1 


-Jj-  =  i^kJ.  -  pk'^  -  v  -  V  ■  Jk, 


(27) 


or,  Fi'nm  equation  24, 


(20) 


ujhere  ue  have  used 


Ck  =  —  =  concentration. 
P 


(20) 


Conservation  of  momentum  can  be  written  as  a  balance  of 
forces  in  the  form 


av  „  -e^  - 


(30) 


Here  Fk  is  the  external  force  per  unit  mass.  We  could  write 
this  Fk  as 


^^kPk  =  Yl^'^  unit  volume), 


(31) 


but,  De  Groot  is  writing  the  equation  for  a  chemical  reaction 
where  mass  is  a  more  important  variable. 


Conservation  of  energy  is  written  as 


d  /I.  .  ,, 

p —  —V  •  V  +  U 

^  dt  V2 


V  [P 
\ 


^-Ffc  ■  vkpk, 


C3E) 


for  Jq  representing  a  heat  flow  Cusually  electrical  in 
origin)  and  P  is  the  pressure  reduced  to  a  scaler  quantity 
from  its  tensor  quality. 

Gibb’s  equation  is  the  final  piece  of  the  puzzle 
necessary  to  connect  all  of  the  important  state  variables. 


dt  di  dt  di 


C33) 


Again,  we  have  the  state  variables  S  ,  JJ  and  P  are  per  unit 


mass  . 


lie  can  take  the  force  balance  equation  Cequation  30), 
multiply  by  the  center  of  mass  velocity  v  ,  and  subtract  the 
equation  for  conservation  of  energy  to  obtain 


=  -PpV-v-  \^.J^  +  Y^Fk-  fk 


C34) 


This  equation,  plus  the  equation  of  conservation  of  mass- 
Cusing  Ck  ),  can  be  substituted  into  Gibb’s  equation  to  yield 


pT—  —  -  V  •  J,j  +  P,i  ■  Jk  P  ^  Pk'^  ■  Jk  -  dc  ^ , 


C35) 


c. 


De  Groat  rearranges  this  equation  to  explicitly  shouj  the 
thermodynamic  driving  forreo  uji<-h  their  conjugate  fluxes 


dS  „ ^  Jq  ■  Jk  ■  Xi, -r  AJ c 

'’a  - T - )' - T - 


=  -V  J 


(.2 


liihere 

X-u,  =  ~^VT  ,  the  thermal  driving  force, 

=■  ,  the  diffusion  force, 

A  =  ,  the  chemical  affinity, 

J,  =•  the  entropy  floui, 

=  entropy  production. 

We  can  now  a'~^ive  at  the  simpler  form 


+  V  •  SvV, 


where  S„  is  the  entropy  per  unit  volume. 
Equation  36  can  be  reformed  as 


dt 


Js  +  S^v]  +  E. 


C3 


C3 


Now,  the  total  entropy  flow  contains  a  convective  term  also 


Equation  30  is  in  the  form  of  a  continuity  equation  for 
entropy  flow.  The  entropy  production  term,  3;  .comes  from 
the  irreversible  processes  that  occur  as  heat  conduction, 
diffusion  or  chemical  reactions.  Since  the  velocity  v  is 
not  in  E  ,  the  bulk  motion  of  the  system  is  considered  a 
reversible  phenomenon. 

We  noui  apply  equations  36  and  30  to  a  system  where 
electrical  conduction  occurs  to  determine  if  the  continuity 
equation  for  entropy  changes  form.  Since  no  chemical 
reaction  occurs  we  can  drop  the  chemical  affinity  from  the 
equations.  The  external  forces  due  to  an  applied  potential 
or  electric  field  in  a  solid  are 


Fk  =  qkE  ^  -ek^<P- 


Here,  the  electric  potential  is  ^  and  et  is  the  charge 
Ccharge  per  unit  mass!  on  the  kth  component  of  concern.  In 
solid  the  applied  potential  also  alters  the  chemical  poten¬ 
tial  or  Fermi  level.  Thus,  the  chemical  potential  takes  the 
form 


The  internal  energy  must  also  include  the  added  electrical 
energy  of  the  carriers  as 


MV.-CV.V-J.N-  V 
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U  -  ^■'  ■  ekCirfp. 


The  thermodynamic  Flux,  or  energy  fiouj  J,  ,  must  now  include 
a  contribution  from  the  current  inhere 


y  ~  ckJk<p  =  14), 


C4e) 


and  I  is  the  total  current  density.  Following  De  Groot’s 
procedure,  just  outlined  to  arrive  at  the  entropy  rate 
equation,  uie  start  with  balance  of  forces 


dv  — -k  -  ■  ^  *— ♦ 

^li  ~  (^kCk^4)  -  V-  P  . 


The  Fk  term  contains  the  external  non-electrical  Forces. 
The  energy  equation  becomes 


d  - 

p —  -V  ■  V  )-  U 

^dt\2 


-V  ■  (P  v-^r  P 


1441 


^Fk  -r-kPk  -  P  l^^kCk-^ 
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Gibb*s  equation,  with  eieccrical  contributions,  becomes 


di  dt  '  dt 


Edc  ^ ^  d(j) 

k  k 


C45) 


Following  De  Sroot’s  procedure  we  arrive  at 


.F'\7  .  ,7  -  V  • 


+  V'Fk  -fk  +  YskCk  —  . 


C46:) 


Substituting  equation  45  and  the  convective  derivative 


for  the  concentration  into  Sibb’s  equation  results  in  the 
following  form 


dS  _  ^  /  Jci  -  UkJk  Jq  ■  Flu  +  Yk  Jk  ■  Fk 
^  dt  ~  \  T  '  T 


C47) 


and 


X,  =  A  -  TV 


\T  J 


=  Xk  -  ek(f>Xu 


B1 


C48) 

Xk  -■  Fk  ~  ek^(f> 


This  is  expanded  further,  using  equation  48,  as 


Jq  Y^k  t^kJk 
\  '  T 


Jq  ■  A\  +  Jfc  •  X;t 

T 


L48) 


which  is  the  form  arrived  at  previously.  The  forces  and 
fluxes  derived  for  electrical  conduction  are  invariant  under 
a  linear  transformation  of  fj.^ ek<P  • 

Analusis  of  Ridleu's  Derivation 

Uie  examine  Ridley’s  derivation  to  insure  that  his 
procedure  did  not  yield  a  result  different  then  De  Groot’s 
results.  Ridley’s  derivation  starts  with  Gibb’s  equation 
where  the  variables  have  been  expressed  as  per  unit  volume. 


y.vy. 


dS  _  db  dpk 

di  dt  dt 

k 


C50) 


For  a  solid,  the  volume  does  not  change  under  current 
conduction,  uihich  allous  us  to  drop  the  work  due  to  pressure 
and  volume.  The  use  of  pk  comes  From  multiplying  De  Groot’s 
previous  form  of  Gibb’s  equation  by  p  to  change  the  variables 
to  per  unit  volume. 

Ridley  introduces  the  additional  electrical  contribu¬ 
tions  of  the  chemical  potential  and  interral  energy  as 


U  =  U  +  ^  ekpk4> 


C51) 


—  >Ik  +  ek<f> 


Which  is  the  same  as  De  Groot’s  formulation.  The  energy 
balance  equation  can  be  uiritten  as 


dU  d 

+  j, 


C52) 


-V  •  (  ^  tkJk^  1  +  ^2 


This  form  matches  De  Groots  formulation  for  the  change  in 
internal  energy  . 


Ridley  substitutes  equations  51  and  58  into  Gibb’s 
equation  to  obtain  a  form  similar  to  De  Groat’s  formulation 


dt 


-V  ■ 


ft} 


T 


Jo 


^  ^  l^kJ k  ‘ 
k 


+f-j  +  Y,x,-j\-Y, 

k  k 


C53) 


ujhere 

F  =  -V</)  , 

J  ”  Sfc  ^kJk  , 

A%  =  -X„  , 

Xfc  =  . 

Ridley  has  broken  the  ohmic  term  F-f  out  explicitly  in 
equation  53,  and  has  allowed  for  local  fluctuation  of  the 
mass  density,  an  aspect  De  Groot  did  not  allow.  The  formula¬ 
tion  still  demonstrates  that  the  change  in  entropy  is  a  sum 
of  the  negative  divergence  of  flow  and  a  positive  production 
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Conditions  oF  the  Stationaru  State 

We  can  nouj  apply  the  conditions  of  the  stationary  state 
to  our  entropy  rate  equations  to  determine  the  action  of  the 
state  variables  in  the  solid. 

Prigogine  C68D  describes  a  stationary  state  as  a  condi¬ 
tion  where  the  state  variables  do  not  depend  on  time.  Yet, 
the  net  mass  flow  is  zero,  and  entropy  production  does  not 
vanish.  This  is  considered  a  non-equilibrium  state  of  an 
open  system  as  opposed  to  an  equilibrium  state  of  an  isolated 
system  where  there  is  no  entropy  production  Cbecause  all 
forces  vanishl . 

In  a  sense,  the  non-equilibrium  stationary  state  is 
similar  to  a  steady  state  of  the  system  with  the  condition 
that  the  entropy  production  rate  is  not  zero. 

Prigogine  CB3J  characterizes  the  stationary  state  as  a 
situation  where  we  allow  the  transfer  of  matter  and  energy 
between  two  phases  of  a  system  at  different  temperatures  such 
that  the  entropy  production  rate  is  positive,  or 


dS 

—  J th^th  Jm^m  ^  0- 


□nsager’s  equations  allow  us  to  view  the  fluxes  and 
forces  in  terms  of  linear  relations  such  that 


C54) 


C55D 


Jm  -  T-2\^lh  't-  L-n^m- 


For  the  stationary  state,  the  mass  flow  is  zero 


Jm  —  L'^i^th  +  J'22^m  =  0. 


C55D 


C5B) 


Qnsager’s  reciprocal  relations  require  L12  =  J-'ti  that 


^  =  LuXf,  +  2L2lXtKXm  +  U2XI  >  0. 
di 


C57) 


If  uje  take  the  derivative  with  respect  to  at  constant 


Xth  we  arrive  at 


d  /dS 


dX^  \  dt 


=  2^,t2lAt/i  +  L22Xmj  —  2J-n„  —  0.  C58) 


Therefore,  the  two  conditions 


T  r  ^  r. 


C53) 


are  equivalent  as  long  as  Unsager’s  equations  are  of  linear 
form  for  the  system. 

The  procedure  followed  above,  where  the  stationary 
states  are  characterized  by  the  conjugate  flows  being  zero. 
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is  similar  to  the  analogous  mechanical  situation  where  the 
conjugate  momentums,  of  the  Forces,  vanish.  The  stationary 
state  defines  a  minimum  for  the  entropy  production  rate  and 
is  referred  to  as  the  state  of  minimum  entropy  production. 

If  a  stationary  state  exists,  then  the  entropy  is  independent 
of  time  and  the  entropy  production  rate  must  be  balanced  by 
the  outward  flow  of  entropy. 


Application  of  the  Stationaru  State 
bu  Ridleu  to  the  Case  of  CCNR 


In  this  section  we  look  at  how  Ridley  applied  the 
concept  of  the  stationary  state  to  his  entropy  rate  equation. 
Ridley  applies  the  idea  of  the  stationary  state  to  the 
situation  of  CCNR  and  VCNR  by  using  the  following  argu- 
mentC701 . 

c  S 

In  the  steady  state  '2i  Cthe  total  entropy  change! 
is  zero,  as  are  the  overall  rates  which  appear  in 
equations  like 


dt 


=  -V  • 


-r 


^kPk 


dcj) 

dt 


and 


dpk 

di 


-V  •  j\  + 


dpk 

dt 


and  in  each  case  the  production  is  balanced  by  the 
divergence  of  the  flow.  At  equilibrium  both  the 
production  and  flow  are  zero.  Since  it  is  to  be 
expected  that  the  steady  state  will  be  as  near  the 


equilibrium  state  as  the  various  constants  uiill 
alioui,  ue  expect  the  ejntropy  production  rate  will 
be  as  small  as  possible.  This  plausible  argument 
is  at  the  basis  of  the  principle  of  least  entropy 
production . 

The  equilibrium  Ridley  talks  about  is  a  thermostatic 
equilibrium;  all  the  driving  Forces  vanish.  This  is  net  the 
case  for  the  situation  he  iiiishes  to  portray.  UJhy  he  has 
stated  this  is  not  clear  and  it  may  have  contributed  to  later 
criticisms  of  the  paper. 

Foliouiing  Prigogine’s  arguments,  .Ridley  assumes  uie  Tiave 
a  stationary  state;  no  mass  floui  and  a  balance  of  entropy 
flouj.  Ignoring  heat  flows,  the  remaining  term  for  the  total 
change  of  entropy  contains  an  entropy  production  rate  which 
should  be  at  a  minimum  value  for  a  stationary  state. 

From  equation  53,  the  remaining  term  for  the  entropy 
production  rate  is 


fJ'kJk  ■  Xq  ^  F  ■  J 


V  _  V  — 

^  n'  ^  Tk  dt  • 


C50) 


The  Tk  term  allows  different  carrier  species  to  have 


different  temperatures.  In  a  stationary  state  no  mass  flow 
occurs  so  Jk  is  zero.  But.  since  the  charge  Flow  is  in 
f i lamentsCf or  CCNR)  and  layersCfor  VCNR3  the  mass  flow  is  not 
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not  zero  for  these  structures.  Furthermore,  Ridley  assumes 
that  there  is  no  temperature  or  chemical  potential  gradients 
in  these  structures.  This  implies  that  the  driving  forces  Xg 
and  A'l,  are  zero.  Since  Ridley  assumes  these  structures  are 
stationary  states,  the  rate  of  change  of  the  structures  at 
their  interfaces  must  be  equal,  but  opposite. 

These  arguments  leave  the  joule  heating  term  as  the  sole 
source  of  entropy  production  in  the  structures  or 


T-E  =  F-f=F-Y,eJk. 

k 


CBl) 


Perturbation  Under  CCNR 

For  the  case  of  CCNR  Ridley  explains  the  curve  in  figure  IS 
as  a  single  filamentary  structure  with  cross-section,  a  , 
which  is  a  percentage  of  the  total  cross-section. 


T £.  E\  •/'i  (1  —  o,)  -{-  E2  J2 O'! 


CBE) 


where 


Figure  19.  Current-voltage  characteristic 
exhibiting  current  contralled  negative  resist¬ 
ance.  The  figure  is  not  drawn  to  scale. 


91 


"t"  6J2O'  —  0. 


CB4) 


Ule  nouj  apply  a  perturbation  through  the  field  variables 
as  required  by  Prigogine’s  results  for  equation  set  59. 


El  =  Eq  -r  ~  •  E2  —  Eo  -r-  8E2 


T  -  T  .  r 

—  Jo  — 

dll, 


J2  —  •.^0  "i"  ■  j-i  8  E2 


(.66) 


to  arrive  at 


a  dE  ^ 


C67) 


When  ^  is  negative,  as  seen  in  experiment,  the 
dE 

production  rate  decreases  and  filamentary  structures  are 
favored.  This  is  in  accord  with  Prigogine's  statement  that 
the  contribution  of  the  rate  of  change  of  forces  to  the  total 
entropy  production  is  negative  or  zero  in  a  stationary  state. 
But,  Prigogine  C64]  states  that  this  is  true  only  under  the 


following  conditions 


llHWM.il 


ifvxv-  vair. 
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1 .  Linear  phenomenaiogical  laws  apply  . 

E.  Onsager’s  reciprocal  relations  are  valid. 

3.  Phenomenological  coefficients  can  be  treated  as 

constants . 

Violation  of  these  conditions  are  the  basis  for  criti¬ 
cism  of  Ridley’s  paper  by  Volkov  and  Kogan  C59] ,  and  Takeyama 
and  Kitahara  CBOD .  Let  us  examine  Ridley’s  arguments  to  see 
if  he  has  violated  Prigogirie’s  conditions. 

An  examination  of  condition  one,  in  light  of  the 
current-voltage  characteristic  of  figure  19  certainly  seems 
to  be  a  flaw  in  Ridley’s  reasoning  based  on  the  fact  that 
CCNR  and  VCNR  are  inherently  non-linear  in  their  entirety. 
But,  examination  of  Ridley’s  work  shows  that  his  analysis 
applies  to  the  negative  differential  region  only.  Experimen¬ 
tal  evidence  at  that  time,  and  later  do  not  dispute  that  this 
region  is  linear  in  a  number  of  cases.  Thus,  Ridley’s  paper 
does  not  fall  into  dispute  due  to  che  first  condition. 

If  the  negative  differential  region  is  linear,  and  no 
extensive  heat  build-up  is  allowed,  then,  quantities  depen¬ 
dent  on  heating  are  not  changing.  If  we  proceed  further  into 
the  positive  resistance  region,  heat  build-up  is  quite 
likely,  from  Ridley’s  development,  and  linearity  is  lost. 
Thus,  Ridley’s  derivation  satisfies  the  second  and  third 
conditions  required  by  Prigogine’s  arguments. 

It  is  interesting  to  observe  that  the  two  papers  cited 
also  try  to  utilize  Prigogirie’s  methods  by  linearizing  the 
particular  phenomena  they  believe  is  causing  negative 
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resistance  in  a  material.  Both  papers  shouj  that  CCNR  occurs, 
but  only  if  the  generation  and  recombination  mechanisms 
present  are  linear  Ca  fact  that  is  not  clear  from  experi¬ 
mental  evidence).  Ridley’s  arguments  are  far  more  general. 

In  his  view,  ujhat  ever  the  mechanism,  it  must  be  linear  in 
the  negative  resistance  regime  or  the  negative  resistance 
regime  cannot  occur  Cit  ujould  violate  the  total  entropy 
production  rate  for  a  stationary  state)  .  Any  proposed 
mechanism  would  have  to  he  demonstrated  to  be  linear  in  order 
to  apply  it  to  the  stationary  state  conditions. 

As  shown  earlier  in  chapter  II,  and  other  studies  C551, 
negative  resistance  can  be  caused  by  a  transition  from  a 
single  carrier  current  regime  to  a  double  carrier  current 
regime.  As  double  carrier  current  begins  to  dominate  current 
flow,  the  negative  resistance  egion  disappears. 

Therefore,  as  the  current  transitions  via  CCNR,  we  can 
treat  the  negative  resistance  regime  as  a  stationary  state. 
This  allows  us  to  utilize  the  general  results  of  Ridley’s 
work  as  we  model  current  flow  under  CCNR  as  a  solid  state 


plasma  . 
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LHAPTER  IV 


QUANTITATIVE  ANALYSIS  OF  THE  SOLID 


STATE  PLASHA  PINCH 


In  chapter  III,  ue  deterrr.ined  that  the  conditions 
surrounding  CCNR  place  limits  of  action  on  certain  macrosco¬ 
pic  variables  dealing  with  driving  forces  that  influence 
carrier  transport. 

If  uie  model  the  solid  state  plasma  under  these  limita¬ 
tions,  it  should  be  possible  to  determine  a  set  of  parameters 
that  describe  the  plasma’s  motion  under  CCNR.  Application  of 
a  perturbation  to  the  plasma’s  motion  under  CCNR  uill  allow 
us  to  determine  the  plasma’s  stability. 

The  motion  of  a  plasma,  under  the  influence  of  electro¬ 
magnetic  fields,  are  best  handled  using  the  equations  of 
motion  for  a  fluid  and  naxweli’s  equations.  Utilizing  this 
information,  attention  in  this  chapter  will  first  be  given  to 
the  similarity  between  diffusive  fluid  flow  and  diffusion  of 
magnetic  fields  in  preparation  to  solve  for  the  equilibrium 
condition  between  the  thermal  and  magnetic  forces  that  occur 
in  and  around  the  solid  state  plasma  under  CCNR  conditions. 

This  same  problem  will  then  be  solved  using  a  pertur¬ 
bation  to  determine  if  the  magnetohydrodynamic  equations  are 
stable  under  CCNR. 
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The  Basic  riLdrodunamic  Equations 


The  basic  moment  equations  for  conservation  of  mass. 


momentum  and  energy  are 


I  V  .  (,t?)  =  0, 


=-^V?+i(/xB), 


Dt  3  cr  Di  ap^ 


where 


P  =  the  mass  density, 


U  =  the  velocity, 


P  ^  the  pressure. 


J  =  the  current  density. 


g  =  the  magnetic  field  strength, 


a  =  the  conductivity. 


7  =  the  ratio  of  specific  heats. 


flaxwell’s  equations  can  be  written  as 


U*J^ 

l^>- 


S' 
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c  c  \ 


Taking  the  curl  of  the  last  equation  yields 


_  _  -  4. TO"  /1_  ,77  t:.  1  dB 

V  X  V  X  5  =  - ( -  V  X  (?7  X  5) - ;— 

;•  V  c  c  ai 


and  using  V  •  B  =  0  we  obtain 


Attct  dB  ^.0-  47rcr_ 

—r~r-  -  ^  ^ - X  {V  X  B). 

at 


By  defining  the  magnetic  viscosity  as 


47rcr  ’ 


we  can  rearrange  equation  71  to  obtain 


^  ^  V  X  ([/  X  5), 

dt 
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The  First  term  on  the  right  hand  side  is  a  diffusion  term  and 


the  second  term  is  a  convection  term.  This  equation  is 
similar  in  form  to  the  hydrodynamic  equation  of  floui  past  an 
object . 

t 

If  the  fluid  is  at  rest,  the  convective  term  vanishes  | 

and  uue  arrive  at  a  solution  similar  to  Stoke ’s  equation  1I66D  > 


dB 

di 
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Equation  74.  like  Stoke ’s  equation,  relates  the  diffusion  of 
physical  quantities.  From  a  fluid  context,  the  diffusion 
equation  states  that  the  rate  of  fluid  density  changes  by  the 
transfer  of  molecules  from  one  region  to  another  Ci.e.,  the 
process  of  diffusion).  For  the  hydromagnetic  case,  the  rate 
of  change  of  the  magnetic  field  changes  by  diffusing  through 
the  electrically  conducting  fluid. 

Continuing  the  similarities  between  fluid  and  hydromag¬ 
netic  diffusion,  uie  can  observe  that  the  magnetic  viscosity 
can  be  related  to  how  fast  rhe  field  diffuses  in  some  time  t 
(similar  to  a  lifetime)  into  the  fluid  or 


AttctL^ 

where  L  is  a  characteristic  length  Cor  penetration  depth) . 
Convection  can  dominate  diffusion  in  the  same  manner  as 


inertial  forces  in  a  fluid  can  dominate  viscous  forces.  A 


measure  of  the  dominance  of  convection  to  diffusion  is  found 


in  the  ratio  of  these  values  that  yields  the  Reynolds  number 
CB71 
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m 
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The  conditions  of  the  stationary  state  Cfrom  chapter 
III)  cause  diffusive  forces  to  vanish.  This  implies  the 
solid  state  plasma  is  highly  conductive  compared  to  the 
semiconductor  surrounding  it,  or,  the  Reynold’s  number  of  the 
solid  state  plasma  is  much  smaller  than  unity  under  CCNR 
conditions.  Ule  can  utilize  this  information  to  solve  for  the 
balance  of  forces  in  the  plasma. 

The  Static  Pinch 

The  similarities  between  fluid  and  hydromagnetic 
descriptions  allow  us  to  look  at  the  effect  of  a  magnetic 
field  on  the  motion  of  a  highly  conductive  fluid.  In 
particular,  we  can  examine  rhe  effect  of  the  pinching  term, 
under  the  conditions  where  the  fluid  is  at  rest  with 
the  magnetic  field  lines.  This  does  not  mean  the  fluid  is 
not  flowing  in  the  material,  merely  that  the  magnetic  field 
is  frozen  into  the  fluid  and  moves  along  with  it  Cas  a  conse¬ 
quence  of  the  fluid’s  high  conductivity). 

The  procedure  will  start  with  the  hydrodynamic  equation 
for  conservation  of  energy  in  inviscid  flow.  Equations 
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relating  the  velocity  of  waves  parallel  and  perpendicular  to 
the  field  lines  in  the  fluid  will  be  used  to  arrive  at  a 
simple  relation  between  the:  outward  thermal  pressure  and  the 
inward  magnetic  pressure.  This  relation  will  be  found  to  be 
the  Bennett  pinch  criterioii. 

Under  the  assumption  of  a  small  Reynold’s  number  we 
have  what  is  called  steady,  creeping  flow.  Under  these 
conditions  the  velocity,  U  ,  is  inviscidC independent  of 
time) .  From  reference  71  tne  hydrodynamic  equation  for 
conservation  of  energy  under  inviscid  flow  is 


V{P  +  ~prj^)  +  p{QxU)  =  0^ 


C77) 


where  the  vorticity  is  defined  as  w  =  V  x  .  Expanding 


the  gradient  of  ^  we  get 


V(^^J  =  t/  X  (V  X  [/  +  ([/  •  V)I7). 


Simplifying  equation  77  yields 
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VP  -t-  p(P  ■  \’)U  =  0 


C79) 


Using  conservation  of  momentum  with  the  rate  of  change 
of  velocity  being  zero  yields 
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p(J7  •  V)^  :=  -VP  +  -  Jx  P  =  -VP  +  i-  ^  ^  ^ 

c  Itt 
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P 


+  P )  +  ■  V)t/ -  =  0 

47r  4-77 


47r 


This  is  simplified  further  by  using  the  follouiing 
quantities 


-^{U  -  K),  CBl) 


p 

(47r/7)^/2 


Here,  U  is  the  fluid  velocity  Cuhich  uill  vanish  under 
inviscid  conditions)  and  velocities  V  and  Ul  are  resultant 
vector  velocities  found  from  adding  the  fluid  velocity  U  and 
some  phase  velocity  of  waves  parallel  to  the  magnetic  field 


lines  . 
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The  last  quantity,  Va  >  in  equation  set  B1 ,  is  known  as 
the  Alfv^n  velocity.  This  quantity  is  similar  to  the  wave 
velocity  along  a  stretched  wire.  For  our  case,  the  Alfven 
velocity  is  the  phase  velocity  of  a  wave  along  a  magnetic 
line  of  force.  Substituting  equation  set  81  into  80  yields 


+  —  j  +  .  V)W  +  {W  •  V)F]  =  0, 


CBE) 


and  for  equal  zero  Cinviscid  flow) 


+  ■^)  ~  pi'^o.  •  V)Va 


0. 
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Equation  82  is  the  magnetohydrodynamic  equation  for  inviscid 
flow  and  has  the  same  form  as  the  hydrodynamic  equation  for 
steady  fluid  flow 


VP  +  p{U  ■V)U  =  0. 


C84) 


Clearly,  the  pressure  in  equation  88  has  been  augmented  by 
the  magnetic  pressure,  and  the  Alfv6n  velocity  is  the 
characteristic  velocity  for  waves  along  the  magnetic  field 
lines  when  static  conditions  prevail. 

For  the  stationary  state,  the  current  is  not  changing 
rapidly.  This  implies  a  magnetostatic  situation  may  exist 
for  times  short  compared  to  the  time  CCNR  occurs.  For  a 
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magnetostatic  situation,  uj-iere  the  magnetic  force  is  balanced 
by  the  pressure  gradient,  and  the  fluid  velocity  is  zero, 
conservation  of  momentum,  equation  60,  allows  us  to  write 


p{U  •  V)C7  =  --VP  + 


(V  X  5)  X  5 


cas: 


Where  the  left  hand  side  of  equation  05  vanishes  as  a 
consequence  of  the  stationiiry  state.  From  our  Stuke's  like 
equation  74  we  have 


cas) 


U)e  can  rewrite  the  second  term  on  the  right  hand  side  of 
equation  05  for  the  static  case  as 


(V  X  j?)  X  J9  1  -  1  _  /BP 

47r  47r  •  gji- 
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And  in  tensor  notation 


d 

/B,Bk\ 

d 
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\  Ttt  / 

dvk  ' 
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The  riaxwell  stress  tensor  is  defined  in  Jackson  CBOl  as 


‘  SB,  -  +  £=)!. 
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Dropping  the  electric  field  stress  term  ue  arrive  at  a  form 
equivalent  to  that  obtained  for  equation  SB  Cue  are  not 
discounting  the  electric  field  at  this  point). 


„  -B.Bt  B\ 
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At  \  _  /(vxP)  xi 

dvk  drk  \  47r  /  dr^  ^  Btt  ’’’^ )  \  Air 


Our  static  equation  89  nou  reads  as 


CBl) 


^-[P^ik  -  -'Uj  =  0  or  ~F,k  =0. 
OTk  drk 


uhere  Fit  is  defined  as  the  total  stress  tensor  equal  to 
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P^Bk 
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This  can  be  reduced  to  diagonal  form  via  transformation 
to  the  principal  axes 


1*5*1' 

m 


X 


This  reduces,  after  some  algebra  to 


A.  =  P 
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A  geometric  interpretation  of  this  matrix  shows  that  -^3 


is  parallel  to  the  field  and  Aj  ,  A2  are  perp 
field  as  shown  in  figure  20. 


endicular  to  the 


Under  static  conditions,  the  self  magnetic  force  takes 


on  the  form 
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Figure  20.  Geometrical  interpretation  of 
equation  36.  A  tube  of  force  in  a  magnetic 
field  is  like  a  string  under  tension.  Increasing 
the  field  strength  stretches  the  tube  of  force, 
and  the  ujave  velocity  increases. 
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or 


B -vP  r=  f.^P  =  0.  C30) 

This  demonstrates  further  that  the  magnetic  field  lines  and 
the  current  lie  on  lines  of  symmetry  as  draun  in  figure  SO. 
For  a  straight  column  these  equations  give  the  appearance  of 
a  shellCor  shells)  of  current  and  their  accompanying  field 
lines  nested  around  the  ceriter  of  the  pinch  column. 

Ue  can  solve  for  the  seif  magnetic  force  of  the  plasma 
column  in  cylindrical  coordinates  to  obtain 


B  d  , 


C39) 


the 


Hultiplying  by  r"  and  integrating  out 
pinch  Cinitially  at  t’o  )  yields 


to  the  radius  of 


r^o 
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{rB)d{rB), 
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I  f  the 
energy 


thermal  pressure  vanishes  at  the  surface,  and  the 
of  holes  and  electrons  are  the  same,  then 


¥ 


2KbT  7i.2rdr  = 

Jo  O'l 


j'Y  =  2  /  TTrndr  =  carriers/  unit  length, 
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2NKsT  = 


Solving  for  the  current  related  to  the  field 


1  d  47r 
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Solving  equation  103  ue  arrive  at  the  Bennett  pinch  criterion 


2NKbT  =  --  or  =  4NKbTc\ 
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This  result  shouis  hauj  the  macroscopic  variables,  such 
as,  the  current,  temperature  and  carrier  density  are  related 
tc  achieve  a  balance  between  thermal  forces  and  magnetic 
forces  under  stationary  state  conditions  in  the  solid  stats 
plasma . 

This  result  is  also  utilized  in  classical  gaseous  plasma 
discharges  for  determining  the  minimum  conditions  necessary 
to  pinch  the  column. 

At  this  point,  objections  may  be  raised  concerning  the 
previous  derivations  which  lead  to  the  Bennett  pinch  criter¬ 
ion.  In  chapters  I  and  II,  it  was  stated  that  semiconductors 
undergoing  second  breakdown  have  very  intense  electric  fields 
and  heating  present.  Also,  it  is  unclear  that  the  situation 
in  the  device  can  be  approximated  by  a  magnetostatic  case 
during  second  breakdown. 

These  criticisms  are  valid  if  one  considers  the  entire 
second  breakdown  current  characteristic.  Yet,  Chapter  III 
demonstrated  that  a  stationary  state  is  linear  in  character 
during  the  CCNR  regime.  Thi.s,  over  times  that  are  short 
compared  to  the  total  time  when  CCNR  occurs,  the  field  does 
approximate  a  steady  state,  and  under  CCNR  the  plasma  column 
will  reach  a  static  pinch. 

The  Hudromagnet ic  Stability  of  the  Pinch 

Under  CCNR  a  solid  state  plasma  column  appears  to  obtain 
an  equilibrium  condition  described  before  as  a  static  pinch. 


lirUinil7WV(rinir\w^rwxrw  4. 


We  DDUJ  relax  the  conditions  on  the  solid  state  plasma  by 
alloujing  all  of  our  pertinent  variables  to  vary  slightly  from 
the  equilibrium  situation.  If  the  static  pinch  is  in  a  state 
of  stable  equilibrium,  a  perturbation  of  the  pinch  uiill 
return  to  its  equilibrium  condition.  If  the  pinch  is  not 
stable  upon  perturbation,  the  pinch  represents  an  unstable 
equilibrium.  This  can  be  determined  by  solving  for  the 
dispersion  relationship  for  the  pinch  column. 

To  determine  if  an  instability  occurs  uie  will  resolve 
the  magnetohydrodynamic  equations  Cequation  set  6B)  under  a 
small  perturbation.  Variables  uith  an  o  subscript  are  first 
order  quantities  and  all  perturbations  are  noted  with  primes. 


P  =  Pc  +  p', 


I7  =  ll  +  V'  .=  t/', 


P=Po-^P\ 
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5  =  s'  +  B'. 
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Ule  uiill  also  maintain  our  result  From  chapter  III  that 
the  plasma  column  has  a  conductivity  much  higher  than  the 
surrounding  semiconductor.  Ignoring  products  of  perturba¬ 
tions  tue  find  that  the  zeroth  order  quantities  for  equation 
set  68  are  related  as 


VPo  =  -(X>  X  B), 
c 


Pop.-  =  constant , 


ClOB) 
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The  first  order  results  are 


dp' 


dU'  .  1 

Po-^  =  -V.P'  +  -[r  X  Bo  +  Jo  X  5'), 
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P’Po  =  yp'Po, 
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Using  equation  sets  107  and  lOB,  uie  want  to  arrive  at  a 
force  equation  involving  some  small  displacement  of  the 
plasma  column  in  terms  of  zeroth  order  quantities  after 
eliminating  the  first  order  results.  Eliminating  p'  and  /' 
from  the  first  order  equations  yields 


C103) 


for  the  first  order  continuity  equation.  The  displacement  of 
the  plasma  can  be  defined  as 


t)  ~  r  -  To, 
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Integrating  the  continuity  equation  in  equation  set  107 


yields 
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d isplacemRnt  and  tirriR.  The  r.nrrnai  method  of  tackliny  this 


equation  is  expanding  ^  as  a  Fourier  series  of  the  form 


f  =  y  icxp{iu}r,t). 
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Equation  113  nouj  appears  as  folloujs  uith  equation  IIH 


substituted  on  the  left  hand  side  to  yield  a  set  of  norma¬ 


lized  equations 


=  Q(0) 
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ujhere  Wn  are  the  normal  frequencies.  For  >  0  all  modes 


should  be  periodic  or  oscillating  around  some  initially 


stable  point.  If  < 0  the  perturbation  grows  exponentially 


and  the  resulting  configuration  is  unstable 


For  our  boundary  conditions,  we  can  use  the  external 


fields  Coutside  the  plasma)  and  the  fields  from  the  perturbed 


quantities  to  arrive  at 
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Luhcrc  prirr.sd  qucnt iticir  ar3  tlie  results  of  perturbations  suol. 
that  V  X  V  X  yl'  vanishes. 

The  tangential  field  uan  be  Found  from  the  image  force 
of  t.he  plasma  at  the  plasma's  surface  as 


C117) 


nx[E[  +  -{U'  +  Bl,)) 


c 


where  quantities  with  an  i  subscript  are  at  or  near  the 
plasma  column  surface  Cas  a  consequence  of  the  plasma 
column’s  conductivity  being  much  higher  than  the  surrounding 
semiconductor  1  .  Ule  can  then  conclude  that 


E[  +  --({/'  X  Bli)  =  0. 


ciia) 


This  allows  us  tc  wrice  equation  117  with  equation  llu 
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Integrating  equation  119 


A.  —  n  X  X  -Boeit)  —  (n  ■  ^)Bogxt- 
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lohich  is  arrived  at  by  having  the  plasma  column  parallel  to 


the  outer  surface  of  the  semiconductor .  This  is  also  a 


consequence  of  the  static  pinch  results.  The  plasma  column 


initially  filled  the  semiconductor  sample,  and  as  a  result  of 


balance  of  forces,  retreated  a  constant  distance  from  the 


sample  surface.  No  isotropies  in  position  are  allouird 


initially  under  stationary  state  conditions. 


Hydrodynamic  considerations  for  conservati  n  of  energy 


require  the  total  pressure  to  be  continuous. 


{Po  +  VP.)  +  P'  +  (f  •  V)p;,  -f  B\)f 
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Using  our  equilibrium  condition 


Bh  _  Bl. 
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equation  ISE  becomes  upon  linearization  Cusing  our  perturbed 
quantities ) 

C124) 

=  %^(S7. +(f-V)S<,„,) 

^  /I 

after  using  P' =  ~7Po'^  ■  (  —  i ■  V Po  •  now  require  further 

information  on  the  magnetic  field  to  proceed  further. 

IniLialiy  there  is  an  azimuthal  magnetic  field  external 
to  the  plasma.  Under  the  conditions  of  the  static  pinch,  the 
magnetic  field  ujill  not  penetrate  the  plasma  to  any  depth  of 
consequence.  Further  perusal  of  the  procedure  used  to  solve 
a  similar  situation  in  classical  gaseous  plasma  columns 
demonstrates  that  an  axial  magnetic  field  is  necessary  for 
stability  of  the  pinch.  Since  some  experimental  results  shoui 
that  an  axial  magnetic  field  is  important  to  stable  some 
solid  state  plasma  instabilities  (see  chapter  II),  ue  will 
add  an  axial  magnetic  field  to  our  derivation  also.  Later, 
this  condition  can  be  removed.  Our  magnetic  field  conditions 
external  and  internal  to  the  plasma  column  are  nouj 

Boe.t  -=  [0,Bo{r),0) 
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^0,1  ^  {'0,0,  Boz  =  constant). 

To  Further  simplify  the  derivation,  ue  will  assume  that 
no  pressure  gradient  exists  radially  in  the  plasma  column,  as 
required  from  the  stationary  state.  This  may  appear  to  be  an 
oversimplification,  uhen  it  seems  quite  probable  that  a 
concentration  gradient  may  exist  under  a  perturbation.  Yet, 
there  is  no  experimental  evidence  describing  what  type  of 
concentration  gradients  exist,  if  any.  Also,  the  stationary 
state  does  not  allou  any  initial  concentration  gradients 
(concentration  gradients  mould  cause  diffusive  forces,  mhich 
vanish  under  the  stationary  state) . 

Placing  these  conditions  into  equation  1E4  results  in 

=  7i^oV(V  •  0  T  ^[V  X  V  X  (f  X  X  bI,. 

For  cylindrical  perturbations  me  have 

^{r']exp{iTn9  -r  ikz) 


{^r{'r),ie',-),f,z{r))exp{ime  +  ikz). 
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and  since  ujb  are  using  very  Icuj  orders  of  approximation  for  j 

our  static  pinch,  we  will  only  utilize  m  equal  zero.  Results 
for  the  classical  plasma  pinch  indicate  m  equal  zer^  corre¬ 
sponds  to  inducing  a  sausage  instability  and  m  equal  one  a 
kink  instability  CBSi  .  Researcn  documented  by  Snyder  C70J 
suggests  that  semiconductor  devices  could  exhibit  plasma 
instabilities  resembling  a  sausage  or  kink  instabilities 
depending  on  the  actual  geometry  of  current  flow.  Expanding 
123  with  m  =  0  in  equation  128  yields 
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Using  the 


definitions  c,  = 


1/2 


P  / 


for  the  speed  of  sound,  and 


1/2 


■i^a  —  l4-p]  for  the  Alfv6p  velocity,  uie  can  simplify  these 


equations  to 
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These  equations  relate  the  displacement  of  the  pinch  to  the 
perturbation  initially  started  with  equation  set  105. 
Eliminating  Cr  from  the  above  equations  to  arrive  at  the 
displacement  in  terms  of  the  z  component  yields,  after  some 
algebra , 
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This  equation  is  in  the  form  of  a  Bessel  equation,  and 
the  solutions  are  modified  Bessel  functions,  (z  =  Io{kr)  .  The 
/„  are  not  solutions  at  r  equals  zero.  Our  argument  for  the 
Bessel  function  is  rather  complicated  and  takes  the  extended 
form 


'W^Tr.'’ 


1  "V*  '-h  .-v>  "J*  “>  V»  "V**  ' 


ri!K‘ 


I 

*^.  *>  *> 


kJclvl  -a;2(c2  +vl) 


.efi  ,  (t)' _ 

clvl  -  (f )'  {c]  +  vl) 


C133) 


and  solving  For  the  z  component  of  the  displacement 
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If  we  use  5'  =  V  X  X  Bo)  for  B'-.  the  z  component  is 
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Placing  the  z  components  into  our  perturbation  result, 
equation  1S4,  ue  get 
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From  n  X  A'  =n  y.  {(  X  Boe)  =  •  i)Boe 


uie  see  that 
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where  a  is  same  arbitrary  function  so  that 
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These  last  two  equations  imply  that 
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—  =  /(7-)  and  -  —{i^Be)  =  g{z). 


But,  since  ivBe  is  a  function  of  r  only  g[z)  must  be 
constant  and 


V  X  A'  =--  g{z)e  =  0, 


since  E'  is  zero  initially.  Ue  can  now  rewrite  equation 
to  get  closer  to  the  dispersion  relation 
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■  This  situation  is  stable  if  an  axial  magnetic  Field  is 
present  such  that 


I-IrKiHr) 


Io{Hr) 


C147) 


We  can  expand  the  dispersion  relation  further  by 
approximating  the  Bessel  functions.  Standard  expansions  of 
these  functions  are 
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If  uje  expand  them  to  first  order  only,  our  dispersion 
relation  Cequation  144)  sirnnlifies  to 
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utilizing  the  definition  for  the  speed  of  sound  in  a  material 


ujhere 
kn  - 


Debye  wavenumber , 


=  plasma  frequency , 

and  eliminating  the  axial  magnetic  field  leaves  our  first 
order  dispersion  equation  as 
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For  this  solution  to  be  stable  mould  require 
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mhicb  is  not  a  physically  admissable  solution  in  the  materi¬ 
al  . 

These  results  shorn  a  static  plasma  pinch  is  unstable 
under  the  conditions  of  CCNR .  The  groming  oscillations  of 
the  pinch  mould  certainly  result  in  anistropies  in  carrier 
concentration  in  the  column,  resulting  in  a  movement  amay 
from  the  stationary  state. 

This  action  alloms  thermal  and  electrical  driving  forces 
to  reappear  and  the  rate  of  entropy  production  mould  rapidly 
increase,  and  possibly  surpass  the  flom  of  entropy  from  the 
sample.  This  suggests  rapid  ohmic  heating  from  an  decrease 
in  the  conductivity  of  the  plasma  column.  Groming  oscilla¬ 
tions  at  the  column  surface  mould  certainly  be  limited  by 
phonon  dispersion  mhich  mould  transfer  the  kinetic  energy  of 
plasma  oscillations  to  the  lattice,  aiding  heating  of  the 
device.  An  exact  determination  of  this  type  of  activity  is 
beyond  the  scope  of  this  monk,  but,  is  consistent  mith 
results  that  shorn  devices  exhibit  signs  of  thermal  shock 


after  the  initiation  of  second  breakdomn . 
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Solid  State  Plasma  Parameters 
Necessaru  to  Initiate  Pinching 

The  results  of  the  last  section  indicate  attempts  to 
determine  when  a  device  uiiil  fail  after  second  breakdown 
occurs  is  an  attempt  to  model  an  inherently  non-linear  event. 

Therefore,  it  seems  more  amenable  to  determine  when  a 
CCNR  regime  occurs  under  device  operation.  The  Bennett  pinch 
criterion  could  be  utilized  to  determine  a  minimum  condition, 
beyond  which,  self-induced  magnetic  forces  overcome  thermal 
forces  and  produce  unstable  pinching.  This  can  be  related  to 
the  quantum  degeneracy  parameter  to  yield  a  set  of  conditions 
that  are  inherently  unstable  under  CCNR  conditions. 

Our  starting  point  is  the  Bennett  pinch  criterion,  as 
applied  to  a  situation  where  equilibrium  is  initiated.  In 
this  case,  the  thermal  forces  Just  balance  the  magnetic 
forces  .  Further  increases  in  the  magnetic  force  initiates 
unstable  pinching.  This  can  be  quantified  as 
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where  the  current  is  in  statamps  and  r  is  the  column  radius. 
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This  value  can  be  substituted  into  our  relations  for  the 
quantum  degeneracy  parameter,  or  the  plasma  coupling  para¬ 
meter,  to  determine  a  set  of  maximum  current  curves  for  a 
material . 

Results  are  shown  for  the  materials  of  interest  on 
figures  El,  EE  and  E3  using  material  parameters  from  EB4D . 
Values  of  N  are  graphed  as  concentrations  per  cubic  centi¬ 
meter  and  all  results  are  for  a  plasma  column  of  radius  10"^ 
centimeters . 

For  the  lowest  current  shown  on  the  figures,  the 
equivalent  current  density  is  greater  than  3  A/cm} 

To  determine  the  current  required  to  support  an  unstable 
plasma  pinch,  the  concentration  of  carriers  in  the  pinch  is 
followed  horizontally  across  the  graph  until  it  intersects 
the  appropriate  isotherm.  These  values  are  placed  into 
equation  155  to  determine  the  maximum  current  that  can  flow 
before  magnetic  forces  destabalize  the  plasma.  For  example, 
on  figure  El,  a  concentration  of  around  6.3  x  at  300K, 

requires  about  10“^  amperes  to  become  unstable  with  a  column 
radius  of  10~^ cm  .  Additionally,  the  quantum  degeneracy 
parameter  at  this  point  indicates  the  plasma  is  classical  in 
nature,  which  reaffirms  our  results  which  indicated  pinch 
results  similar  to  those  found  in  classical  gaseous  plasmas. 

These  results  also  indicate  that  for  a  given  concentra- 
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tion,  temperature  and  current,  more  than  one  filament  could 
form.  The  smaller  filaments  would  each  have  to  posses  values 
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Figure  22.  Graph  of  the  concentration 
versus  the  quantum  degeneracy  parameter  for 
pinch  conditions  in  silicon.  Pinch  radius  is 
2Q-2  centimeters. 
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Figure  23.  Graph  of  the  concentration 
versus  the  quantum  degeneracy  parameter  for 
pinch  conditions  in  gallium  arsenide.  Pinch 
radius  is  centimeters. 
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belauj  the  ma.'iinum  allcuiable  values  fer  the  Bennett  pinch 
criterion,  and  could  be  treated  as  a  set  of  current  carrying 


ujires  with  appropriate  repulsive  forces.  As  the  current  in 
the  material  increases,  filaments  ujould  coalesce  back  into  a 
single  filament  uihich  mould  become  unstable  once  it  achieves 
the  maximum  allowable  current  from  the  Bennett  pinch  criter- 


Further  results  are  indicated  on  figures  22  and  23  for 
silicon  and  gallium  arsenide  respectively. 
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CHAPTER  V 


APPLICATION  OF  THE  PLASIIA  PARAflETERS  TO 
SEN  I  CONDUCTOR  JUNCTIONS 

In  this  chapter,  the  Bennett  pinch  criterion,  derived 
for  a  solid  state  plasma  in  a  CCNR  regime,  will  be  applied  to 
a  transistor  that  is  placed  into  second  breakdown  and 
recovered  safely  before  damage  can  occur.  These  results  for 
the  transistor  are  the  only  known  experimentally  repeatable 
results  of  second  breakdown  and  provide  a  means  to  determine 
if  second  breakdown  in  this  device  coincides  with  the 
conditions  for  unstable  solid  state  plasma  pinch. 

As  an  introduction  to  these  results,  a  short  discussion 
on  the  current  transport  mechanisms  in  solid  state  devices 
Ce.g.,  diodes  and  bipolar  transistors)  will  be  presented. 

The  results  of  this  discussion  will  be  used  to  determine 
whether  or  not  the  transistor  to  be  examined  in  a  CCNR 
regime.  Measurements  of  the  collector  current  at  second 
breakdown  will  be  used  to  determine  the  carrier  concentration 
available  for  a  plasma  whose  initial  radius  is  equal  to  the 


When  two  dissimilar  semiconductors  are  brought  into 
contact,  a  built-in  potential  is  set  up  due  to  the  difference 
in  Fermi  energies.  This  built-in  potential  sweeps  free 
carriers  from  the  general  area  of  contact,  in  both  materials, 
until  enough  charge  has  been  displaced  to  counteract  the 
built-in  potential.  The  remaining  charge  free  zone  is  called 
a  depletion  region. 

Application  of  an  external  potential  changes  the 
electrical  character  of  the  depletion  region  depending  on 
whether  the  potential  aids  Ccalled  reverse-bias)  or  counter¬ 
acts  Ccalled  forward-bias)  the  depletion  region’s  built-in 
potential . 

A  single  Junction  between  two  semiconductors  is  usually 
called  a  diode.  In  reverse-bias  a  diode  passes  very  little 
current  until  the  external  potential  is  large  enough  to 
initiate  impact  ionization  in  the  depletion  region.  As 
reverse-bias  increases,  the  electrical  character  of  the 
depletion  region  changes,  from  an  insulator,  to  a  semiconduc¬ 
tor,  when  impact  ionization  occurs  Csee  chapter  II).  In 
forward-bias  the  built-in  potential  is  decreased  by  the 
external  potential  allowing  a  large  number  of  carriers  to  be 
injected  and  a  large  current  to  flow.  Under  forward-bias 
conditions  the  depletion  region  is  small  enough  that  current 
flow  resembles  semiconductor  injection  Csee  chapter  II). 


If  three  dissimilar  semiconductors  are  brought  into 
contact,  with  the  middle  semiconductor  relatively  thin,  tu/o 
depletion  regions  occur  separated  by  a  very  small  distance. 
Injection  activity  at  either  Junction  can  influence  injection 
conditions  at  the  other  Junction,  creating  uihat  is  called, 
the  transistor  effect. 

Under  different  polarities  and  degrees  of  bias,  these 
Junctions  can  be  used  as  an  electronic  switch,  also  known  as 
a  transistor  Csee  figure  24).  When  the  transistor  is  on, 
both  Junctions  are  forward-biased  Cthe  transistor  is  in 
quasi-saturation)  or  the  base-emitter  is  forward-biased  and 
the  collector-base  is  slightly  reverse-biased  Cthe  transistor 
is  active).  By  reverse-biasing  both  junctions,  neither 
Junction  will  pass  current  and  the  transistor  is  off  Cin  the 
cut-off  state).  Both  situations  are  portrayed  in  figures  24 
and  25. 

These  effects  are  attributed  to  the  control  of  charge  in 
the  base  region  by  the  emitter’s  bias,  and  the  degree  of  base 
drive  used  to  drain  or  replenish  stored  charge  in  the  base. 

Under  active,  quasi-saturated  or  cut-off  conditions,  the 
quantity  of  stored  charge  varies.  As  shown  in  figure  26,  the 
ability  of  a  transistor  to  switch  Cfrom  an  on  to  an  off 
stateJ  IS  controlled  by  how  quickly  charge  is  drained  out  of 
the  base  Cor  the  collector)  to  achieve  cut-off  conditions. 
This  is  complicated  if  the  external  bias  on  the  collector  or 
emitter  continues  to  pour  charge  into  the  base  region. 


Figure  24.  Schematic  depiction  of  a  CNPN) 
transistor.  The  upper  figure  shows  the  extent 
Cnot  to  scale)  of  the  depletion  regions  with  the 
built-in  potential  only.  The  bottom  figure 
shows  the  transistor  Cnot  to  scale)  in  a  quasi- 
saturated  state.  Both  depletion  regions  are 
reduced  by  forward-bias. 
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depletion  regions  (cut — off  state) 


Figure  26,  ninority  carrier  density  stored 
in  the  various  regions  of  a  NPN  transistor  in 
various  conduction  states.  Here,  n^,  is  the 
electron  density  injected  into  the  p  type  base. 
The  hole  density  being  injected  i’~*o  the  n  type 
collector  is  p,-,  .  Very  little  charge  is  stored 
in  the  emitter. 
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For  example,  from  the  results  of  chapter  II,  a  transis¬ 
tor  could  be  modeled  as  a  diode  and  PIN  diode  connected  in 
series.  Special  consideration  must  be  taken  into  account  in 
modeling  the  intervening  base  region,  but  current  flows 
outside  the  base  will  still  meet  the  same  primary  obstacle  of 
space  charge  limited  current  flow. 

Conditions  can  be  created  where  the  collector  and 
emitter  potentials  force  their  respective  depletion  regions 
to  extend  deep  into  the  base  until  they  coalesce  Cknown  as  a 
punch-through  condition)  .  Once  this  occurs  a  transition  from 
space  charge  limited,  single  carrier  to  double  carrier 
current  flow  occurs  and  continued  current  flow  cannot  be 
restricted . 

If  conditions  are  created  during  transistor  turn  off 
where  carriers  are  depleted  from  the  base  at  the  expense  of 
carriers  in  the  collector,  space  charge  limited  current  flow 
in  the  collector  can  eventually  increase  the  electric  field 
at  the  I-N+  interface  Csee  figure  25)  to  avalanche  condi¬ 
tions,  Avalanche  creates  a  large  number  of  electron-hole 
pairs  which  neutralize  the  space  charge  and  leads  to  double 
carrier  current  flow  and  CCNR  C71)  . 

These  situations  provide  a  useful  test  of  the  conclus¬ 
ions  reached  in  chapter  IV.  Knowing  the  current  levels 
available  at  the  collector  prior  to  second  breakdown,  the 
collector  area,  the  charge  carrier  concentration,  and  the 
maximum  thermal  velocity  of  carriers  in  the  collector  region 


uie  can  see  if  the  minimum  predicted  conditions  for  solid 
state  plasma  pinch  coincide  with  experimental  evidence  of 
second  breakdown  initiation. 

The  only  repeatable  experiments  driving  transistors  in 
to  second  breakdown,  without  damage,  have  been  performed  by 
Portnoy  C7S,73I1.  In  the  next  section  we  will  discuss  these 
results  and  apply  our  pinch  criterion  to  the  tested  devices. 

Application  of  Plasma  Pinch  Criterion 
to  Transistors 

Portnoy’s  experimental  set-up  examines  transistors  in 
the  common  emitter  configuration  as  they  are  switching  under 
a  large  inductive  load  C723 .  When  second  breakdown  occurs, 
the  inductive  load  is  quickly  switched,  from  the  device,  to 
an  energy  dump,  preventing  device  failure.  We  will  use  the 
results  of  testing,  listed  in  reference  73,  for  a  transistor 
manufactured  by  Unitrode,  with  general  characteristics 
described  in  figure  27. 

Results  of  this  study  demonstrate  that  three  types  of 
second  breakdown  occurred  dependent  on  the  device’s  initial 
condition  and  the  degree  of  reverse  base  drive.  These  types 
of  second  breakdown  have  been  called  types  A,  B,  and  C  C7E1 
by  Portnoy . 

Utilizing  information  derived  in  chapter  IV,  a  graph  of 


the  necessary  pinch  conditions  for  the  device  were  calculated 
for  a  column  radius  of  0.2B  centimeters  Csee  figure  20). 


Figure  E7 .  Doping  profile  for  the  Unitrode 
power  transistor.  The  area  of  the  device  is  a 
square  0.56  centimeters  on  each  side.  From 
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The  constant  current  lines  on  the  graph  oF  Figure  E8 
represent  maximum  allowable  values  oF  current  at  speciFic 
values  oF  the  temperature  and  carrier  concentration. 

We  can  study  each  type  oF  second  breakdown  individually 
to  determine  iF  CCNR  characteristics  are  being  initiated  and 
compare  the  current  and  charge  carrier  concentration  to  those 
predicted  For  pinch. 

We  cannot  necessarily  speciFy  whether  a  particular 
breakdown  was  initiated  by  punch-through  or  avalanche 
breakdown  at  the  collector.  Both  events  occur  under  very 
similar  conditions  and  represent  the  only  mechanisms  possible 
in  a  transistor  at  these  levels  oF  voltage  and  current. 

Since  inFormation  after  breakdown  is  not  available  Cwithout 
damaging  the  devices)  deductive  evidence  For  either  situation 
is  not  available.  In  either  case,  the  carriers  transiting 
the  collector  will  be  at  their  maximum  thermal  velocity, 
which,  with  the  collector  current  and  area,  allows  us  to 
specify  a  minimum  carrier  concentration  in  the  collector. 


Type  A  Second  Breakdown 

The  general  current  and  voltage  trends  with  time  for 
type  A  second  breakdown  at  room  temperature  are  shown  in 
figure  2S . 

The  initial  state  of  the  transistor  in  type  A  breakdown 
is  active  or  quasi-saturation.  This  implies  the  emitter  is 
forward-biased  and  the  collector  is  forward-biased  or 
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Figure  29.  Typical  waveforms  of  current 
and  voltage  for  type  A  second  breakdown .  The 
top  curve  is  the  voltage  across  the  collector- 
emitter  ClOO  volts  per  square),  the  middle  curve 
is  the  collector  current  C2  amperes  per  square), 
and  the  lower  curve  is  the  base-emitter  voltage 
CO. 5  volts  per  square).  The  time  base  is  5 
microseconds  per  square.  The  waveform  cut-off 
at  the  far  right  is  where  second  breakdown 
occurs.  From  C72i . 
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slightly  in  reverse-biased.  In  either  case,  the  depletion 
regions  are  not  very  large  and  do  not  extend  across  the 
entire  base  region  initially.  The  bias  of  the  Junctions 
creates  an  excess  of  carriers  in  the  base  from  the  forward- 
biased  emitter  Cand  collector).  If  the  collector  is  slightly 
reverse-biased,  a  small  number  of  these  excess  carriers  will 
be  drained  off  by  the  collector. 

As  the  base  drive  is  reversed  to  turn  the  transistor 
off,  the  inductor  senses  the  bias  change  on  the  transistor 
and  begins  to  de-energizes  across  the  transistor  forcing  the 
collector  into  reverse-bias  and  trying  to  keep  the  emitter 
forward-biased . 

General  trends  indicate  that  for  type  A  breakdown  the 
collector-emitter  voltage  remains  constant  and  rises  Just 
before  second  breakdown.  The  collector  current  decreases 
linearly  to  zero  and  presumably  becomes  negative,  and  the 
base-emitter  voltage  decreases  more  or  less  linearly,  then 
drops  rapidly  near  second  breakdown  Csee  figure  29) . 

In  terms  of  the  device’s  switching  dynamics,  the  emitter 
cannot  become  reverse-biased  until  all  the  excess  charge  in 
the  base  region  is  removed  by  the  base  drive  and  the  collec¬ 
tor.  If  this  is  correct,  then  second  breakdown  would  not 
occur  until  all  the  base  charge  is  removed.  Once  emptied, 
punch-through  or  avalanche  at  the  I-N+  interface  can  occur. 


Calculations  of  the  minimum  carrier  concentration  in  the 
collector  are  based  on  a  maximum  thermal  velocity  of 
Results  are  listed  in  table  1 . 

The  table  lists  the  initial  device  conditions,  in  terms 
of  the  forward  base  drive,  and  the  conditions  as  the  transis¬ 
tor  is  being  turned  off  by  the  reverse  base  drive.  The  time 
to  second  breakdown  is  the  time  after  initiation  of  the 
reverse  base  drive  until  a  rapidly  increasing  current  flow  is 
detected,  and  current  is  shunted  away  from  the  transistor. 
This  is  the  transition  from  a  high  voltage,  low  current 
state,  to  a  low  voltage,  high  current  state,  or  second 
breakdown . 

The  tabulated  results  show  the  collector  current  exceeds 
the  maximum  pinch  criterion  for  the  minimum  concentration  of 
carriers  in  the  collector. 

Tupe  B  Second  Breakdown 

The  state  of  the  transistor  for  this  type  of  second 
breakdown  is  quasi-saturation.  Typical  current  and  voltage 
curves  with  time  are  shown  in  figure  30. 

Figure  30  shows  collector  current  is  constant  after 
applying  reverse  base  drive  with  an  initially  positive 
emitter-base  voltage  and  a  small  collector-emitter  voltage. 
Near  second  breakdown  the  emitter  base  voltage  tends  to  zero 
and  the  collector-emitter  voltage  rises  rapidly  to  a  sustain¬ 


ing  yoltage. 
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Figure  30.  Typical  uaveForms  of  current 
and  voltage  For  type  B  second  breakdown.  The 
top  curve  is  the  collector  current  CE  amperes 
per  square),  the  middle  curve  is  the  base- 
emitter  voltage  CO. 5  volts  per  square),  and  the 
bottom  curve  is  the  collector-emitter  voltage 
ClOO  volts  per  square).  The  time  base  is  500 
nanoseconds  per  square.  Waveform  cut-off  at  the 
far  right  is  where  second  breakdown  occurs. 

From  C7E3 . 
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The  initial  conditions  for  type  B  second  breakdown  are 
quite  different  from  type  A.  With  both  junctions  well  into 
forward-bias  the  base  and  collector  regions  are  saturated 
with  accumulated  charge.  As  the  inductor  tries  to  drive  the 
collector  to  reverse-bias,  charge  drains  out  of  the  collec¬ 
tor.  The  emitter-base’s  forward  bias  pours  additional  charge 
into  the  base.  This  minority  charge  injected  into  the  base 
may  have  enough  energy  upon  entering  the  base  to  transit  the 
narrow  depletion  region  and  aid  in  diminishing  the  collec¬ 
tor’s  stored  charge. 

As  the  collector’s  stared  charge  diminishes,  the 
collector  moves  to  reverse-bias.  Since,  the  base  drive  also 
removes  charge,  once  the  collector  becomes  reverse-biased  the 
emitter  also  becomes  reverse-biased  and  the  conditions 
quickly  move  to  punch-through  or  avalanche  as  described  in 
type  A  breakdown. 

Table  B  indicates  that  conditions  are  above  the  maximum 
levels  dictated  for  pinching. 


Type  C  Second  Breakdown 

The  initial  state  of  the  device  is  the  same  as  type  B 
second  breakdown.  Typical  current  and  voltage  curves  with 
time  are  shown  in  figure  31.  The  major  difference  between 
type  B  and  C  second  breakdown  is  the  occurrence  of  a  negative 
voltage  for  the  base-emitter  before  the  collector  became 
reverse-biased . 
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Type  6  second  breakdoujn . 


Figure  31.  Typical  waveform  of  current  and 
voltage  for  type  C  second  breakdown.  The  top 
curve  is  the  collector  current  C2  amperes  per 
square),  the  middle  curve  is  the  base-emitter 
voltage  CO. 5  volts  per  square),  and  the  lowest 
curve  is  the  collector-emitter  voltage  ClOO 
volts  per  square) .  the  time  base  is  500 
nanoseconds  per  square.  Waveform  cut-off  at  the 
far  right  is  where  second  breakdown  occurs. 

From  C7E1 . 
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If  the  emitter-base  becomes  reverse-biased  due  to  a  high 
reverse  base  drive  then  the  collector’s  stored  charge  Cuihich 
is  much  smaller  than  the  base’s  stored  charge)  will  quickly 
be  dissipated  allauing  the  collector  to  move  into  reverse- 
bias  and  punch-through  or  avalanche  conditions  more  quickly 
than  in  type  B  breakdoiun . 

Table  3  shows  once  again,  the  collector  current  far 
exceeds  the  maximum  allowable  pinching  current. 

From  the  quantitative  point  of  view  taken  for  these 
results,  the  criterion  for  solid  state  plasma  pinch  appear  to 
be  satisfied  in  all  cases,  and  the  devices  present  character¬ 
istics  that  indicate  the  presence  of  a  CCNR  regime. 
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All  currents  and  limes  ate  Irorn  {tJO].  Maximum  allowable  pinch  currents 
are  loimd  using  the  cor’centration  in  ths  collector  as  a  minimum  value. 


Table  3,  continued 


CHAPTER  VI 


GENERAL  CONCLUSIONS  AND  FUTURE  WORK 

This  ujork  has  established  the  phenomena  seen  during 
second  breakdouin  Cin  particular,  current  f i lamentation)  in  a 
semiconductor  device  is  a  consequence  of  the  solid  state 
plasma  pinching  under  CCNR  conditions. 

To  accomplish  pinch,  certain  requirements  in  the 
semiconductor  must  be  met .  These  requirements  can  be 
generally  determined  through  the  use  of  the  carrier  concen¬ 
tration,  quantum  degeneracy  parameter,  current,  and  tempera¬ 
ture  as  dictated  by  the  Bennett  pinch  criterion. 

The  solid  state  plasma  pinch  was  calculated  to  be  an 
unstable  event,  under  first  order  approximation,  which  leads 
to  growing  oscillations  and  non-linear  behavior. 


Future  U)ork 

The  calculations  for  stability  of  the  static  pinch  could 
be  solved  to  higher  order  modes  with  different  distributions 
of  carrier  temperatures.  This  is  not  a  simple  task,  and  at 
some  point,  approximation,  or  advanced  numerical  techniques, 
will  have  to  be  utilized. 
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Yet,  it  seems  inevitable  that  this  task  could  be  greatly 
simplified  if  repeatable  experimental  information  on  the 
solid  state  plasma  pinch  were  available. 

This  may  be  accomplished  by  utilizing  microuiave  inter¬ 
ferometry  on  a  semiconductor  sample  undergoing  plasma  pinch. 

This  technique  mould  utilize  comparison  of  the  phase 
shifted  micromave  signal  passing  through  the  solid  state 
plasma,  to  an  unchanged  reference  signal,  to  find  the 
plasma’s  index  of  refraction.  This  information  mill  detail 
the  concentration  of  the  plasma  in  real  time  as  the  plasma 
grams . 


Applications  to  Device  Failure 

The  information  developed  in  this  effort  mas  initially 
prompted  by  the  desire  to  predict  or  explain  device  failure. 

Certainly,  the  information  detailed  here-in  has  shomn 
testing  devices  to  thermal  failure  cannot  accurately  predict 
individual  device  failure.  A  semiconductor  device  is 
electrically  unstable  before  thermal  failure  and  going  beyond 
this  initial  instability  is  extremely  difficult  to  model. 

Application  of  the  plasma  parameters  developed  in  this 
effort  could  possibly  determine  the  minimum  requirements  for 
plasma  pinch  in  some  devices,  but,  some  device  designs  mould 
initiate  a  kink  instability  before  a  sausage  instability. 
Thus,  further  mork  is  needed  to  determine  a  Bennett-like 
criterion  for  the  kink  instability  in  a  semiconductor. 
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Further  application  to  test  cases  ujould  then  determine 


unstable  operating  areas  for  semiconductor  devices  based  on 
the  base  material  parameters,  temperature,  and  current 
levels.  If  successful,  large  scale  testing  of  devices  would 
not  be  necessary . 

Also,  some  researchers  have  found  that  some  plasma 
pinches  in  semiconductors  can  be  stabilized  by  exciting 
higher  order  perturbations  Cm='3,4,...)  in  the  pinch.  This 
concept  will  have  to  be  explored  for  semiconductor  devices. 

If  found  to  be  feasible,  semiconductor  devices  could  be 
designed  to  operate  at  much  higher  current  levels  then  is 
presently  possible. 


LIST  OF  REFERENCES 


1.  Tauc,  S.  and  Abraham,  A.,  Phys .  Rev.  lOB .  pp .  S36-B37 
C1B57)  . 

E.  Thornton,  C.  G.  and  Simmons,  C.  D.,  IRE  Trans.  Elect. 
Dev.  ED-5 .  pp .  6-10  ClSSB)  . 

3.  Snyder,  fl.  E.,  AFUiL-TR-83-SE ,  Air  Force  Weapons 
Laboratory,  Kirtiand  Air  Force  Base,  New  Mexico.  pp .  E- 
3  C13a3). 

4.  Bergman,  F.  and  Geistner,  D.,  IEEE  Trans.  Elect.  Dev. 
ED-13  pp.  630-634  C1365) . 

5.  Kalab,  B.  II .  ,  HDL-TR-1333,  Harry  Diamond  Laboratories, 

U .  S.  Army  Electronics  Research  and  Development  Comand, 
Adelphi,  Maryland,  p.  3  C130E3. 

6.  Schafft,  K.  A.  and  French,  J.  C.  ,  IEEE  Trans.  Elect. 
Dev.  ED-13.  613  C136B). 


7,  English,  A.  C.,  Solid  State  Elect.  6,  511  C1363) . 


8. 

Wunsch,  D. 
NS-15.  QD. 

C .  and 
E44-E53 

Bell.  R.  R., 
(1366)  . 

IEEE 

Trans . 

Nuc , 

.  Sci  . 

3. 

Tasca,  D.  M 

.,  IEEE 

Trans.  Nuc. 

Sci  . 

N5-17. 

pp. 

364-37E 

(1368)  . 

10.  Kalab,  B.  M.,  HDL-TR-1915,  Harry  Diamond  Laboratories, 

U.  S.  Army  Electronics  Research  and  Development  Command, 
Adelphi,  Maryland.  (1380). 

11.  Ward,  A.  L.,  IEEE  Trans.  Nuc .  Sci .  N5-E3 .  pp .  1637- 
1684  (1376). 

IE.  Ward,  A.  L.,  HDL-TR-137a,  Harry  Diamond  Labortories, 

U.  S.  Army  Electronics  Research  and  Development  Command, 
Adelphi,  Maryland.  (138S) . 

13.  Alexander,  D.  R.,  Karaskiewlcz ,  R.  J.  and  Eniow,  E.  W., 
AFWL-TR-8O-1E0  pt  1,  Air  Force  Weapons  Laboratory, 
Kirtiand  Air  Force  Base,  New  Mexico.  (1381). 


170 


Lv  n  L’v  i;v  u’v  uv  'w’wwww 


Morgan,  E.  L.,  Adams,  □.  E.,  Rouian,  Ul .  H.  and  Templar, 

L.  C.  Electrical  Qverstress/Electrostat ic  Discharge 
Symposium  Proceedings.  E05-3 .  pp .  151-166  C19B1) . 

Portnoy,  W.  h .  and  Bamble,  F.  R.,  IEEE  Trans.  Elect. 

Dev.  ED-11 .  470  C1S64:>  . 

Sunshine,  R.  A.  and  Lampert,  f1.  A.,  IEEE  Trans  Elect. 
Dev.  ED-IS.  G7E  C137EJ . 

Snyder,  M.  E.,  AFU/L-TR-85-1 15 ,  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  New  hexico.  C1S87D . 

Schafft,  H.  A.  and  French,  J.  C.,  IEEE  Tran.  Elect.  Dev. 
ED-13,  pp.  613-618  C13B6). 

Lampert,  fl.  A.  and  Mark,  P.,  Current  In  lection  in  Solids 
CAcedemic  Press,  New  York,  New  York,  1970). 

Snyder,  M.  E.,  AFWL-TR-85-1 15 ,  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  New  flexico.  C19B7). 

Jonscher,  A.  K.,  Brit  Jrnl .  Appl .  Phys .  1^,  366 
C1961)  . 

Narsh,  J.  W.,  Mayer,  J.  W.  and  Baron,  R.,  Appl.  Phys. 
Lett.  5,  74  C1965). 

Mayer,  J.  W.,  Baron,  R.  and  Marsh,  □.  J.,  Phys.  Rev. 

137.  aE86  C1365). 

Snyder,  M.  E.,  AFWL-TR -85-1 15 ,  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  New  Mexico.  C1987) 
p.  44. 

Snyder,  M.  E.,  IBID  p.  5E . 

Ridley,  B.  K.,  Proc .  Phys.  Soc .  8E,  954  C1963). 

Alekseev,  M.  E.,  Sov .  Phys.  Semicond .  3,  604  C1970) . 

Osipov,  V.  V.  and  Kholodnov,  Y.  A.,  Sov.  Phys.  Semicond. 
7,  604  C1973). 

Vladimirov,  V.  V.,  Golovinskii,  P.  M.  and  Gorshkov,  V. 
N.,  Sov.  Phys.  Semicond.  3,  1514  C1981). 

Ridley,  B.  K.  and  Prate,  R.  6.,  Phys.  Lett.  4,  300 
C1964)  . 


31.  Langmuir,  I.,  The  Gollected  Uorks  oF  Irving  Langmuir, 
volume  4,  Electrical  Discharge,  volume  5,  Plasma 
Oscillations .  Cpergamon  Press,  Neui  York,  Neui  York, 
C19B1)  . 

3E .  Bennett,  U .  H.,  Phys.  Rev.  4^,  090  C1934). 

33.  Tonks,  L.,  Phys.  Rev.  56,  360  C1958) . 

34.  Eckersley,  T.,  Nature  135 .  104  C1935)  . 

35.  Barkhausen,  H.,  Phys.  Zs .  401  C1919)  . 

36.  Alfven,  H.,  Ark.  flat.  Astro,  and  Phys.  BE9 .  E  C194E)  . 

37.  Rutheman,  G.,  Ann.  Phys.  E,  113  C194B). 

30.  Lang,  UJ .  ,  Optik  3,  333  (1940). 

39.  Pines,  D.  and  Bohm,  D.,  Phys.  Rev.  0^  6E5  C1951)  and 

05,  330  C195E3. 

40.  Steele,  f1 .  C.  and  Glicksman,  h .  ,  Jrnl ,  Phys.  Chem .  Sol. 
0,  E4E  C1959). 

41.  Chynoueth,  A.  6.  and  Hurray,  A.  A.,  Phys.  Rev.  1E3 .  515 
C1961) . 

4E .  Ivanov,  Y.  L.  and  Ryvkin,  S.  fl.,  Sov .  Phys.  Tech.  Phys. 
3,  7EE  C1950)  . 

43.  Galt,  J.  K.  and  Yager,  U) .  A.,  Phys.  Rev.  114 .  1396 
C1959)  . 

44.  Glicksman,  n.,  Phys.  rev.  1E4 .  1655  C1961). 

45.  Buchsbaum,  S.  J.  and  Galt,  J.  K.,  Phys.  Fluids  4, 

1514  C1961)  . 

46.  Glicksman,  n.  and  Steele,  H.  C.,  Phys.  Rev.  Lett.  E, 

461  C1959D. 

47.  Vladimirov,  V.  V.,  Sov.  Phys.  Usp .  1^,  pp .  690-711 
C1975)  . 

40.  Ando,  K.  and  Glicksman,  (1.,  Phys.  Rev.  154  316 
C19671  . 

49.  Gunn,  B.,  Solid  State  Comm.  1.,  80  C19B3)  . 

50.  Vladimirov,  V.  V.,  Sov.  Phys.  Usp.  1J0,  691. 


173 


Chen ,  F .  F . ,  Introduction  to  Plasma  Phusics  and 
Controlled  Fusion,  volume  1.  CPlenum  Press  Nem  York,  Neuj 
York,  19B5)  p.  897. 


Chen,  F.  F. ,  IBID  pp .  1-4. 


Tanenbaum,  B.  9.,  Plasma  Phusics.  ChcGrauj-Hill  book 
Comp.  Neuj  York,  Neu  York,  19G7)  pp .  3-5. 


Ancker- Johnson ,  B.,  Semiconductors  and  Semimetals, 
volume  6,  Inlection  Phenomena,  chapter  11  C1970) . 


LiboFf,  R.  L.,  Jrnl .  Appl .  Phys.  SB,  E530  C1904) . 
Ichimaru,  S.,  Rev.  Hod.  Phys.  1017  C1908). 

Liboff,  R.  L.,  Jrnl.  Appl.  Phys.  SB,  8533  C1904) . 


Kubat,  n.  Pouer  Semiconductors,  chapter  8.  CSpringer- 
Verlag,  Neuj  York,  Neu  York,  1904) 


Volkov,  A.  F.  and  Kogan,  Sh .  h.,  Sov .  Phys.  JETP  85 . 
1095  C19G7) . 


Takeyama,  K.  and  Kitahara,  K.  Jrnl.  Phys.  Sac.  of  Jpn. 
39  183  C1975) . 


DeGroot,  S.  R.  Thermodunamics  oF  Irreversible  Processes 
CNorth-Hol land  Pub.  Neu  York,  Neu  York,  1951)  p.  75. 


Prigogine,  I.  Introduction  to  Thermodunamics  oF 
Irreversible  Processes  CJohn-Wiley  and  Sons,  Neu  York. 
Neu  York,  1951)  p .  94 . 


B3  . 

Prigogine, 

I  . 

IBID,  pp.  75-77. 

G4. 

Prigogine , 

I  . 

IBID,  p.  96. 

G5. 

Snyder,  M. 

E  . 

and  Volkin,  H.  AFWL-TR-05-1 IB .  Air 

Force 

Weapons  Laboratory,  Kirtland  Air  Force  Base,  Neu 
Mexico,  C 1907) . 


Brodkey ,  R.  S.  The  Phenomena  oF  Fluid  Motion  CAddison- 
Wesley  Pub.  Comp.,  Reading,  Mass.  1957)  pp .  3'^-39 . 


Brodkey,  R.  S.,  IBID,  p.  105. 


Jackson,  J.  D.,  Classical  Electrodunamics  CJohn-Wesley 
and  Sons,  Neu  York,  Neu  York,  19B8)  p.  859. 


Boyd,  T.  J.  and  Sanderson,  J.  J.,  Plasma  Dunamics 
CBarnes  and  Noble  Inc.,  Neu  York,  Neu  York,  19G9)  p.  05. 


sWW.Vt 


Snyder,  n.  E.  and  VolFin,  H.  AFWL-TR-BS-l 15 .  Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  Neuj  flexico 
C13B71  pp.  31-37. 

Snyder,  n.  E.,  AFWL-TR-85-115,  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  Neu  Mexico,  C19B71, 
pp.  77-B5. 

Portnoy,  W.  M.,  AFWL-TR-SE-13S ,  Air  Force  Weapons 
Labortory ,  Kirtland  Air  Force  Base,  Netu  Mexico,  C1SB3). 

Portnoy,  W.  M.  and  McMullen,  S.  A.,  AFWL-TR-B4-111 , 

Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 
Neuj  Mexico,  C1SB5)  . 


S' 


(jS: 

I 


tiUUUlllMil. 


